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Показано, что широкий круг явлений в ферромагнитных переходных 
металлах и их сплавах, изучаемый с помощью рассеяния медленных нейтро¬ 
нов, можно последовательно описать с единой точки зрения, используя 
модель Хаббарда и ее обобщения. Обсуждаются некоторые новые возможно¬ 
сти метода нейтронной спектроскопии, которые открываются с вводом мощ¬ 
ных импульсных источников, таких, как ИБР-2. 

И із зЬолѵп Ійаі а питЬег оі ргоЫетз іп Геггошадпоііс ІгапзШоп піеіаіз 
апй іЬеіг аііоуз, зіисііесі Ьу ІЬе пеиігоп зсаНегіпя, сап Ье йезсгіЬес! іп ипіііей 
таппег изііщ ІЬе НиЬЬагй тосіеі апй Да яепегаіігаііопз. А роззіЬіШу о! 
теазигіпк іЬе Зіопег ехсііаііопз іп піскеі изіпд іЬе ЬідЬ Них риізей геасіог 
ІВВ-2 із сіізсиззей. Весаизе о і іЬе ро\ѵег апсі зресіПс сЬагасІегізІісз о! ІВК-2 
іЬе ЬідЬ іпіепзіѵе і1о\ѵ о! пеиігоп’з сап Ье оЬіаіпесІ апсі ІЬе зикйеаіей ехрегі- 
тепі сап Ье регіогтесі іог іЬе геаіізііс Ііте оі теазигешепіз. 

ВВЕДЕНИЕ 

Пятьдесят лет тому назад авторы наиболее полной по тем вре¬ 
менам монографии о магнетизме следующим образом характери¬ 
зовали положение дел в физике магнитных явлений: «Еще недавно 
вопросы магнетизма представлялись исключительно неблагодар¬ 
ной областью для теоретических исследований. Такое положение 
стояло в связи с тем обстоятельством, что внимание исследователей 
направлялось главнейшим образом в сторону ферромагнитных 
явлений, ибо эти последние играли и играют весьма крупную 
роль в технике. Теоретическое же толкование ферромагнетизма 
представляет столь значительные трудности, что и поныне назван¬ 
ная область является одной нз наиболее темных во всей физике» [1]. 

За истекшие десятилетия физика магнитных явлений превра¬ 
тилась в очень обширную и разветвленную отрасль современной 
физической науки [2]. В значительной мере прогресс в изучении 
структурных и динамических свойств магнитных веществ обязан 
достижениям магнитной нейтронографии ІЗ, 4]. Уникальные воз- 
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можности метода рассеяния тепловых нейтронов позволяют полу¬ 
чать информацию о магнитной и кристаллической структуре, рас¬ 
пределении магнитных моментов, о спектре магнитных возбужде¬ 
ний, критических флуктуациях и т. п. Для интерпретации полу¬ 
ченных данных необходимо учитывать электрон-электронные и 
электрон-ядерные взаимодействия в системе и принцип Паули. 
Таким образом, магнетизм можно последовательно описать только 
в рамках квантовой статистической теории конденсированного 
состояния. 

Благодаря наличию внутренней недостроенной пй- или и/-обо¬ 
лочки все свободные атомы переходных элементов сильномагнит¬ 
ные в силу правила Хунда. При образовании кристалла оболочки 
атомов перестраиваются и для ясного понимания свойств кристал¬ 
ла необходимо знать волновые функции и энергии бывших валент¬ 
ных электронов. Расчет энергетических уровней электронов в кри¬ 
сталле весьма сложен [5, 6]; говоря качественно, речь идет о том, 
насколько изменяются при образовании кристалла атомные волно¬ 
вые функции валентных электронов, насколько они делокализу¬ 
ются. 

В теории магнетизма оказался очень эффективным метод, 
модельных гамильтонианов. Без преувеличения можно сказать, 
что огромные успехи физики магнитных явлений в значительной 
мере обусловлены использованием для «теоретического толкова¬ 
ния ферромагнетизма» нескольких очень простых и схематических 
модельных представлений. Этих моделей, которые в настоящее 
время используются в квантовой теории магнетизма, по суще¬ 
ству, три. 

Исторически первой была модель локализованных спинов 
Гейзенберга (и ее предельный случай модель Изинга [7]), затем 
были предложены модель ферромагнитной ферми-жидкости [2] 
и модель Хаббарда [8]. 

Модель Гейзенберга [2—4, 9] основана на предположении о том, 
что волновые функции магнитоактивных электронов в кристалле 
мало отличаются от атомных орбиталей. Считается, что эта модель 
в основном применима к веществам, у которых энергия основного 
состояния отделена щелью конечной ширины от энергии возбуж¬ 
денных токовых состояний, т. е. к полупроводникам и изоляторам. 
Эта модель также хорошо применима к ряду магнитных редко¬ 
земельных элементов, поскольку недостроенные /-слои имеют 
малые эффективные радиусы (для многих редкоземельных /-метал¬ 
лов важную роль в установлении магнитного порядка играет 
косвенный обмен через электроны проводимости). 

Модель коллективизированных электронов [1, 10] или модель 
ферромагнитной ферми-жидкости применяется к металлам, в кото¬ 
рых система бывших валентных электронов образует по всему 
объему кристалла подвижную ферми-жидкость. (Иногда для 
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^-металлов, у которых недостроенные й-слои имеют не столь 
малые эффективные радиусы как /-слои, и потому имеется замет¬ 
ное перекрытие между ближайшими соседями, говорят о смеси 
двух ферми-жидкостей из 8- и гі-электронов.) 

Модель узких зон или модель Хаббарда [8, 11] в определенном 
смысле является промежуточной моделью и первоначально была 
предложена для описания Зй-металлов. В настоящее время модель 
Хаббарда и ее обобщения применяются для описания магнитных 
свойств чистых Зй-металлов и их сплавов, халькогенидов пере¬ 
ходных и редкоземельных металлов, актиноидов и др. 

Таким образом, все три модели по-разному отвечают на вопрос 
о том, как и насколько изменяются волновые функции бывших 
валентных электронов в кристалле. Использование упомянутых 
моделей (или их комбинаций) позволяет описать очень широкий 
круг явлений и получать качественные, а часто и количественные, 
правильные результаты. Иногда (не всегда) очень сложные и тру¬ 
доемкие расчеты электронной зонной структуры практически не 
добавляют ничего нового по сравнению с результатами, полу¬ 
ченными на основе указанных схематических и грубых мо¬ 
делей. 

Однако несмотря на достигнутые успехи в понимании физики 
магнитных явлений из-за трудностей точного расчета электронной 
зонной структуры и проведения всех необходимых экспериментов, 
мы все еще, за исключением редкоземельных элементов, не можем 
с полной определенностью утверждать, какая из указанных мик¬ 
роскопических моделей (или их комбинаций) наиболее адекватно 
приближает реальную ситуацию в том или ином веществе. Вот 
почему определение истинного механизма возникновения магнито¬ 
упорядоченного состояния считается в настоящее время пробле¬ 
мой номер один в теории магнетизма. 

Точное понимание физических процессов, приводящих к появ¬ 
лению ферромагнетизма, особенно важно, когда переходим к тео¬ 
ретическому описанию сплавов ферромагнитных переходных ме¬ 
таллов. Сплавы магнитных металлов широко используются в тех¬ 
нике, поэтому их теоретическое изучение имеет большое приклад¬ 
ное значение. 

В настоящем обзоре покажем, что исследование спектра магнит¬ 
ных возбуждений чистых переходных металлов и их сплавов пред¬ 
ставляет большой интерес для уточнения наших модельных теоре¬ 
тических представлений о природе магнитного состояния в этих 
веществах. Наиболее прямым и удобным методом эксперименталь¬ 
ного изучения спектра магнитных возбуждений является метод 
неупругого рассеяния тепловых нейтронов. В связи с этим обсуж¬ 
даются некоторые новые возможности метода нейтронной спектро¬ 
скопии, которые открываются с вводом мощных импульсных 
источников нейтронов, таких, как ИБР-2. В частности, рассма- 
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тривается возможность прямого измерения стонеровских возбуж¬ 
дений в ферромагнитных переходных металлах. Прямое измерение 
стонеровских возбуждений в широком интервале переданных 
импульсов и энергий имеет большое значение для суждения о сте¬ 
пени локализации или делокализации магнитоактивных электро¬ 
нов в ферромагнитных переходных металлах. Показано, что 
использование модели Хаббарда позволяет с единой точки зрения 
описать весьма широкий круг вопросов магнитного поведения 
чистых ферромагнитных переходных металлов и их сплавов. 


1. МОДЕЛЬ ХАББАРДА 

В последние годы довольно большое распространение при опи¬ 
сании магнитных и электрических свойств переходных металлов 
и их сплавов получила модель Хаббарда [8, 11, 12]. В настоящей 
работе покажем, что спектр возбуждений 
модели Хаббарда и некоторых ее простых 
модификаций представляет значительный 
интерес с точки зрения применения метода 
рассеяния нейтронов. Сначала будем говорить 
о чистых Зй-металлах, а более конкретно — 
о никеле, железе и кобальте. Считается [13], 
что исходный энергетический спектр системы 
представляет собой широкую $р-зону, в кото¬ 
рую погружена система пяти узких пересе¬ 
кающихся сі - зон (рис. 1). В модели Хаббарда 
сделана попытка учесть весь комплекс не¬ 
обычных свойств Зй-металлов. В частности, 
из эксперимента известно, что спин-волновое 
рассеяние медленных нейтронов 13, 4] в этих 
веществах можно описывать на основе мо¬ 
дели Гейзенберга. Найдено также, что рас¬ 
пределения зарядовой плотности в ферромаг¬ 
нитных металлах близки к атомным распределениям [2—4]. С дру¬ 
гой стороны, средние атомные магнитные моменты заметно отли¬ 
чаются от атомных и являются дробными. О значительном вкладе 
й-электронов в низкотемпературную теплоемкость с эл = я 3 /3^(е / )х 
ХТ = уТ говорит то, что коэффициент у у переходных металлов 
больше, чем у нормальных. О сильной коллективизации й-элект- 
ронов говорят также оценки энергии связи и исследование 
ферми-поверхностей. 

В модели Хаббарда [8, 11, 12] эти аспекты поведения системы 
удается описать, предполагая, что й-электроны образуют зону, 
но испытывают сильное кулоновское отталкивание на одном узле 
решетки. 



О дШ 


Рис. 1. Схема зонной 
структуры ферромаг¬ 
нитного переходного 
металла 


7-01068 



37п 


I. Л. КУЗЕМСКИЙ 


Рассмотрим систему й-электронов, гамильтониан которой имеет 
обычный вид [4] 

я =2[+гі+”<ЧІ+ 1/2 2 |7+тг= л + г <‘> 

» іФ; 

Перейдем к представлению вторичного квантования, вводя 
полевые операторы ^ а (г) и *Р$(г). В этом представлении гамиль¬ 
тониан (1) приобретает следующий вид: 

Я = 2 5 (РгЧІ (г) кЧ'о (г) + 

+ 1/2 2 5 ] (ІЬг&г'Ч'І (г е ) ѵ;» (г) 74+ (г') Ч'а (г). (2) 

Собственные функции невозмущенного гамильтониана к образуют 
полную ортонормируемую систему функций, называемых блохов- 
скими функциями , {ср кпа (г)}. Индекс п характеризует какую-то 
данную зону. Используя разложение полевых операторов через 
операторы а кпа и а кпа для блоховских состояний ц> кпо '. 

(г) = 2 ф*па (г) аіпс, (3) 

перепишем гамильтониан (2) в виде 

# = 2 е п (к )а% пя а кп п + 

кпо 

+ 1/2 2 3 З^пт, рН к - к '> Ч)«а поак‘+ Ч 'ро’а к+Ч тоа к 'іо- (4) 

кк'ятп а о' 

Р/ 

Здесь 

Ѵ П т, рИ к , к '> Ч) = <ФЬпо (О ФА+9РО' (Оі^(г—ОІФй+вто'фА'/о)- ( 5 ) 
Собственные функции невозмущенного галшльтониана ц> кпс (г) 
можно разложить в^ряд Фурье следующим образом: 

ФапДг) = N- ^/2 3 Фпо (Г - В,) ехр (ікК+ ( 6 ) 

Функции ф па (г — К г ), известные как функции Ванье, локализо¬ 
ваны в пространстве и ведут себя подобно атомным волновым 
функциям; они образуют полную ортонормированную систему 
функций, полученную из фь па (г) с помощью унитарного преоб¬ 
разования, коэффициенты которого равны ехр (ікК г ). Число ато¬ 
мов в кристалле обозначено N. Точный характер локализации 
функций Ванье проанализировать довольно трудно [5]; ясно, 
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однако, что на больших расстояниях функции Ванье малы. В реаль¬ 
ных расчетах предполагается, что функции Блоха, соответствую¬ 
щие волновому вектору к в п -й зоне, можно вычислить одним из 
известных методов, например, с использованием Ха-метода [6]. 

Запишем теперь гамильтониан (2) в представлении вторичного 
квантования в базисе Ванье: 


-#=2 2 а іто а т'в + 

И тт'о 

+ 1/22 2 2 (і та ,\І п °' | у-1 кпъ’о, Іп'о'У а1 та а1 па 'аі п .„га кт 'о; (7) 

і Ш тп ' со' 


*Г'= 5 ^ Г Фт (г— КД А (г) фт’ (г —К;); 

^ іто , ]па’ | -і- 1 кт'а, Іп'а ^ = 

* (г — к г ) ф* (г'— К г ) ф т , (г—В 6 ) ф (г' —К е ) і 3 гй 3 г' 


(8) 

(9> 


Матричный элемент 1™ т ' описывает перескоки электрона с атома 
на атом для разных орбитальных состояний; чем больше частота 
перехода, тем сильнее коллективизация электронов, т. е. ширина 
полос тем меньше, чем прочнее связаны электроны. 

Для простоты, в этом разделе ограничимся случаем одной 
зоны. Удобно также представить гамильтониан (7) в следующей 
форме: 


Н = 2 ^/ 40+0 + ^72 2 п іа п і-а-\- 

+ 1/2 2 2 (*і | Т | п * п і°' — 1/2 2 2 1 іі а іо а іо'а1о'ауо + 

а аа' и о0' 

+ 2 2 (“|т| Іі) '<40'«+- (10) 


Здесь введены обозначения: 

1] = е г | \ф (г —КД| 2 |г — г'Г 1 |^(г' — К г )| 2 й 3 гй 3 г'; 

1 и = е 2 \ф*( г- К,) ф (г- КД | г - г' I" 1 ф * (г' - КД ф (г' - КД ѴгсРг’ , 

где — интеграл прямого обмена; V — интеграл внутриатом¬ 
ной кулоновской корреляции электронов с противоположными 
спинами. Третий член в (10) с точки зрения магнетизма не пред¬ 
ставляет интереса, так как он пропорционален А 13 п ь п Последний 
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член в (10) описывает переходы антипараллельных спиновых пар; 
эти эффекты не будем рассматривать здесь. Оценки показывают [4], 
что для многих переходных металлов и их соединений наибольшим 
интегралом в (10) оказывается Я. Исходя из этого, Хаббард [8] 
предложил следующую модель: 

(И) 

Отталкивание электронов на одном центре Я не убывает при раз¬ 
ведении атомов. 

Таким образом, в модели Хаббарда заключена возможность 
перехода от простой зонной схемы к описанию типа Гайтлера — 
Лондона. Поэтому можно ожидать, что наиболее интересные эффек¬ 
ты: переход от обычного металла при | | II, описываемого 
зонной схемой, к диэлектрику с сильной корреляцией гайтлер- 
лондоновского типа при | іц | <§; II, а также различные магнит¬ 
ные свойства должны содержаться уже в гамильтониане (11). 
Возможность построения интерполяционных решений — от атом¬ 
ного до зонного предела, описывающих схлопывание щели в спек¬ 
тре одночастичных состояний,— была показана в [14—17]. 

В случае одного электрона на атом (п = 1) при сильном взаимо¬ 
действии (II | Ц] [) в рамках теории возмущений второго поряд¬ 
ка можно показать [11], что эффективный гамильтониан Я э ф = 
= Я (II) (Е о — Я (іі]))' 1 Н (II) сводится к антиферромагнитному 
гамильтониану Гейзенберга: 

іі 

Таким образом, основное состояние системы для сильной корре¬ 
ляции и наполовину заполненной зоны скорее является антиферро- 
магнитным. 

Гамильтониан Хаббарда обладает свойством ротационной инва¬ 
риантности. Это означает, что при действии оператора спинового 

вращения Я = ехр { — ^ іф 3 °А } (Я 1- /? = ЯВ + = 1), где ф — 

угол вращения вокруг единичной оси к справедливо равен¬ 
ство Я' = ВНВ+. Здесь о іх — (я,| ац ■+• а г \ ); а іу = 
= — і (я4 а Н — а і\ а і\)'і — ( п і\ — п н)- Относительно обра¬ 

щения времени Т (Т + Т = ТТ* = 1) можно также показать инва¬ 
риантность гамильтониана Хаббарда: Я' = ТНТ + . Обращение 
времени приводит для кристалла к выполнению важного равенства 
е (к) == е (— к) вне зависимости от пространственной симметрии 
системы. Свойство ротационной инвариантности гамильтониана 
также приводит к важным физическим следствиям общего харак¬ 
тера, к обсуждению которых и перейдем. 
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2. ОБОБЩЕННАЯ ВОСПРИИМЧИВОСТЬ И СПИНОВЫЕ ВОЛНЫ 

Согласно теории линейного отклика [18], можно ввести обоб¬ 
щенную спиновую восприимчивость X (ч, со). Если внешнее маг¬ 
нитное поле достаточно мало, намагниченность, вызываемая им, 
определяется по формуле 

М(т, Л 3 г'Лі'Х{ г-г', I'). (12) 

После фурье-преобразования получим М (ч, со) = X (ч, со) X 
X Я (ч, со). Здесь X (ч, со) описывает отклик системы на внешнее 
магнитное поле Я (г, і) = Н (ч, со) ехр [і (чг — соД]. Если поле 
реально, намагниченность 

М (г, Д = Я(ч, со)Ке{Х(ч, со)ехр[і(чг — соО)}- (13) 
Частотная зависимость восприимчивости определяется запазды¬ 
вающим коммутатором 

Х(г —г', со) = -і- ^ йі ехр (ісоО {\М (г, I), М (г, 0)]). (14) 

Для свободных электронов статическая спиновая восприимчивость 
(восприимчивость Паули): 

Х р = 2рв-0 (еД, 

где р в = еКІ(2тпс) — магнетон Бора; Я (еД — плотность состоя¬ 
ний на уровне Ферми. 

Обобщенная спиновая восприимчивость для системы взаимо¬ 
действующих электронов определяется следующим выражением: 

Х ар (г, г', (о) = (^ в ) 2 ~ | ^ехр(шД<[Я“(г, і), 5 р (г', 0)]). (15) 

Для гамильтониана Хаббарда (И) в приближении Хартри — Фока 
п п п п -► (га и ) щ ; тіц {пц ) получим для энергий электро¬ 
нов в магнитном поле: 

= г к + Я (іц) — (*ц в /2) Я; 

Е к \ = + Я (пД -)- (^р в /2) Я. 

Статическая восприимчивость (15) при Т = 0 в этом случае вычис¬ 
ляется элементарно 

Х = (*»|4/2)Я(вД/[1 -ЯЯ(еД]. (16) 

Этот результат был впервые получен Стонером. Величина 8 = 
= (1 — ЯЯ (еД) -1 получила название стпонеровского усиления. 
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Критерий возникновения магнетизма Стонера имеет вид ЯЯ (е у ) > 
> 1 и связан с возникновением особенности восприимчивости. 
Возникновение особенности восприимчивости (15) при <?, со -ѵО 
в самом общем случае было проанализировано Эдвардсом и Фише¬ 
ром [19]. Поскольку для модели Хаббарда полный спин является 
интегралом движения ([2 Я] =0), то операторы 8 ± = 
= 8" ± і 8ѵ изменяют общую 2 -компоненту на ±1. При этом 
с гамильтонианом коммутирует также оператор Эдвардс 

и Фишер показали, что из свойства ротационной инвариантности 
модели Хаббарда следует наличие спин-волнового полюса у обоб¬ 
щенной спиновой поперечной восприимчивости Х + ~ (ц, со) следую¬ 
щего вида: 

Шд = Пд ,2 + 0(д' 4 ); 

°Я 2 = ТШг{Ья№чі —ЙѴ'ІітІітХ^. (17) 

> &Ы-0 д-0 

Здесь введены обозначения: 

Х/ = Т 1 *вхр(і©*)<[/5(і), 

%/, = [5“(ц), Я]. 

В пределе со -ѵО, д/со ->0, выражение (17) является точной фор¬ 
мулой, справедливой для любого металлического или неметалли¬ 
ческого ферромагнетика или неферромагнетика в статическом 
магнитном поле. В общем случае Г) определяется зонной структу¬ 
рой кристалла. Для простой кубической решетки приближение 
хаотических фаз дает для модели Хаббарда следующий резуль¬ 
тат [4]: 

# = "6Д- 2 №*} (%* + п йД + Ж 2 "^” 1 2 {^Ь 6 *} 2 ( П к\~~ п к\)- (19) 

к к 

Здесь е й =2 ^ехр{ — ік(К* — КД} — зонная энергия; Д = 
= ( 2 ( й рі — ге р Д) — зонное расщепление. Обобщенная воспри- 

р 

имчивость Х + (д, со) удовлетворяет важному правилу сумм: 

| Іш Х + " (д, со) гісо = -у (га* — «Д = — -у- ($*>. 



( 20 ) 
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Точная формула Эдвардса и Фишера [19] для восприимчивости 

х “-тг 1 +-&-{ х '-тГ< № < 2 ‘> 

показывает, что при д = 0 имеем только первый член в (21), кото¬ 
рый отвечает полюсу спин-волнового типа. Правило сумм (20), 
таким образом, выполняется. Ясно, что для малых д вклад спин- 
волнового полюса должен преобладать. Это обстоятельство имеет 
очень большое значение при анализе рассеяния медленных нейтро¬ 
нов. Таким образом, наличие спин-волнового полюса в спектре 
магнитных возбуждений системы является следствием ротацион¬ 
ной инвариантности гамильтониана. 

Еще в ранней работе Херринг и Киттель [20] показали, что 
в простых приближениях спиновые волны равно хорошо можно 
описать как в рамках модели локализованных спинов, так и в- 
рамках модели коллективизированных электронов [21]. Поэтому 
изучение, например, температурной зависимости среднего момен¬ 
та в N1 и Ре в рамках низкотемпературной спин-волновой теории, 
как правило, не дает указаний в пользу той или иной модели. Для 
проверки зонной теории магнетизма предлагали разные способы 
(см., например, [22]). Интересный результат был получен в работе 
Фонера и др. [21] с помощью эффекта Мёссбауэра. Основываясь 
на предсказании Вольфарта [23], они исследовали поведение Ее 
и N1 в очень сильных магнитных полях (выше технического насыще¬ 
ния) при Т = 4,2 К. В данной области вклад спиновых волн 
практически полностью подавлен. Они нашли, что восприимчи¬ 
вость Хнг = дМІдН > 0, в то время как для модели локализо¬ 
ванных спинов при тех же условиях должно быть Хне = 0. Это 
важный аргумент в пользу зонной модели магнетизма. Однако для 
интерпретации эксперимента в [21] используется простейшая 
модель Стонера (16), и в результате выводы носят лишь качествен¬ 
ный характер. 

С нашей точки зрения [24, 25] наиболее отчетливо различие 
между моделями проявляется в спектре магнитных возбуждений. 
Коллективизированная модель имеет более сложный спектр, чем 
модель локализованных спинов. Этот спектр представляет собой 
наибольший интерес с точки зрения применения метода рассеяния 
медленных нейтронов, который уникален при непосредственном 
изучении магнитной динамики магнитоупорядоченных веществ. 

3. ТЕОРИЯ РАССЕЯНИЯ НЕЙТРОНОВ В ПЕРЕХОДНЫХ МЕТАЛЛАХ 

Среди экспериментальных методов, используемых в физике 
твердого тела, нейтронные методы занимают совершенно особое 
место. Рассеяние медленных нейтронов позволяет определить 
кристаллическую структуру, положение атомных магнитных мо- 





376 А. Л. КУЗЕМСКИЙ 


ментов, распределение неспаренных электронов, частотные спект¬ 
ры и кривые дисперсии различных элементарных возбуждений. 
В данном разделе рассмотрим рассеяние медленных нейтронов 
магнитными возбуждениями переходного металла. 

Как хорошо известно [3, 4], сечение неупругого рассеяния ней¬ 
тронов выражается через мнимую часть обобщенной спиновой вос¬ 
приимчивости Х“Р (д, со). Поперечная спиновая восприимчивость 
определяется следующим образом: 


Х + “(д, со) = (ёЩв) 2 | <Я + (д, *) <$"( — Ч)»ехр (Ш)сВ, (22) 
где 

<А (*), Я» = — І- Ѳ (І) {[А (і), В]) 

— двухвременная запаздывающая температурная функция Грина 
[18]. Спиновые операторы в (22) выражаются операторами вторич¬ 
ного квантования в ^-представлении [4]: 

5+ (д) = Р (д) 2 8 + к (я) = Р (д) 3 а + к+ч , 

5- (д) = Р (я) 3 8\ (?) = Р (д) 2 а + кі а к+ч , 

*Чч)= 5 й 3 г ехр (іцг)І^ (г)| 2 . 

Тогда сечение неупругого рассеяния нейтронов на поперечных 
спиновых компонентах запишется в виде [4] 


г {ІтХ*-(Ч. ш)-ІтХ»-(-ч, -»)}. (23) 


Поскольку структура обобщенной спиновой восприимчивости 
и вид ее полюсов определяются выбором модельного гамильтониана 
системы и характером сделанных при ее вычислении приближе¬ 
ний, то по результатам нейтронных экспериментов можно судить 
об адекватности микроскопических моделей. Однако, чтобы надеж¬ 
но судить о применимости той или иной модели, необходимо изме¬ 
рить восприимчивость во всех точках фурье-пространства и для 
всех температур, что не всегда позволяет имеющаяся эксперимен¬ 
тальная техника. Расчет спиновой восприимчивости Х +_ (д, со) 
для модели Хаббарда (11) в приближении хаотических фаз [4, 26] 
приводит к следующему результату: 

ѵ+- /„ _ X- (д, со) _ 

Л (д, со)- 1 _ і//(№в)2ХГ( д, со) • 


( 24 ) 
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ХГ <№,)’ 2 ; (25) 

Пко = <«ка) = {ехр [р ( Е (ко) — еу)] + !}-'; Е (ко) = е к + -^2 п р°- 


Формула (24) является динамическим обобщением результата 
Стонера (16) для Т Ф 0. 

Полюса полной восприимчивости (24) задаются уравнением 


Из уравнения (26) при д, со ->( 
для спин-волнового полюса 



Рис. 2. Возбуждения в коллективи¬ 
зированной модели. 
Возбуждения типа 1 происходят без пере¬ 
ворота спина, а возбуждения типа 2 — 
с переворотом 


), | д |/со ->-0 следует выражение 
= вида (19). Среди полюсов 



0 Ч 1 


Рис. 3. Спектр спин-волновых и сто- 
неровских возбуждений в однозон¬ 
ной модели Хаббарда 


восприимчивости (24) содержатся полюса хартри-фоковской вос¬ 
приимчивости (25): 

ІтХ+-(я, ю) = я(^р в ) 2 ІѴ" 1 2(п и+ < І> * — п кі )6(Па — е к+ч + е к — Д). 

(27) 

Возбуждения, определяемые законом дисперсии 

ймц = е к+ д — е к + Д, (28) 

принято называть стонервескими возбуждениями. В отличие от 
спин-волновых возбуждений переворот спина определяется здесь 
не коллективным движением системы, а одночастичным (рис. 2). 
Схематически вид спектра возбуждений, определяемого уравне¬ 
нием (26), представлен на рис. 3. 
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Вдали от точки пересечения стонеровского и спин-волнового 
спектра можно приближенно разделить вклады в сечение рассея¬ 
ния от этих двух типов возбуждений. Для этого представим 
Іт X (д, со) в виде 

Іт Х + -(д, (о) = ІтХ 0 + -{[1-С7/(^р в ) 2 КеХ;-]2 + 

ЧДС/Д^ХПТ 1 - (29) 

Так как ІшХ+~(д, ю) = 0, когда (д, со) не принадлежит области, 
в которой заданы стонеровские возбуждения, то из (29) следует, 
что 

(Л) 2 1тХ+ -( 4 ’ ю)= — я6{1—Н/(#р в ) 2 КеХ;-(д, со)}. (30) 
Поскольку 

1 -*7/(^ в ) 2 Не Х 0 + - (д, со) (Ш (0) + На) -±-, 
то из (30) получим 

Іт Х + ~ (д, со) = — л (^р в ) 2 б (Йсо + Ш ((})), (31) 


где ±(} = д — т; т — волновой вектор обратной решетки. С уче¬ 
том (29) — (31) получим из (23) сечение спин-волнового рассеяния 
медленных нейтронов: 


й 2 а 

<т іе 


к (2я) 3 Г уе 2 \ 2 
к 0 ѵ 0 ( т е с 2 ) 




х 2 2 (я (О (<») 6 (»® + «2 ГС)) б ГС +д-т) + 

О т 

+ [п (О ГС)) + 1] 8 (На - Ш (0)) б (д- (*- т)>, (32) 

= п (й) = [ехр (РШ)— I] -1 . 


Сечение рассеяния (32) в точности совпадает с сечением рассеяния 
на спиновых волнах для модели Гейзенберга, вычисленным в при¬ 
ближении хаотических фаз (см., например, [4]). 

Однако для модели Хаббарда (11), как это следует из уравне¬ 
ний (27), (29), в области энергий, где электрон в состоянии преодо¬ 
леть энергетический барьер, связанный с переворотом спина 
в эффективном поле, должно наблюдаться рассеяние нейтронов 
на стонеровских модах [3]: 

(-шг) № = (т^ГІт 1*49)141 + Й- 0 >> X 

х 2 2 Яко(1 — гак+ Ч а') б (Я(к-{-д 0 ') — Я(к(т — Па)). (33) 
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Таким образом, неупругое рассеяние нейтронов происходит 
в данном случае вследствие возбуждения переходов | ко > -> 
-► | к + цр' >• Ферромагнитные пере¬ 


ходные металлы имеют хорошо выра¬ 
женные пики неупругого рассеяния на 
спиновых волнах. Для многих метал¬ 
лов закон дисперсии Бф измерен с 
большой точностью. На рис. 4 при¬ 
ведены результаты нейтронных изме¬ 
рений для N1. По своему характеру 
сечение рассеяния на стонеровских воз¬ 
буждениях (33) сильно отличается от 
сечения рассеяния на спиновых вол¬ 
нах. Стонеровское рассеяние согласно 
(33) должно приводить к сильно раз¬ 
мытым диффузным пикам малой интен¬ 
сивности. Поскольку такие возбужде¬ 
ния не возникают в модели Гейзенберга, 
то их непосредственное обнаружение 
с помощью неупругого рассеяния ней¬ 
тронов являлось бы важнейшим аргу¬ 
ментом в пользу модели Хаббарда. Это 



Рис. 4. Закон дисперсии 
спиновых волн в N1 [99] 


тем более актуально, что существует 

большое число косвенных данных о их существовании. Мы об¬ 


судим эти данные в следующих разделах. 


4. МАГНИТНЫЕ СВОЙСТВА НИКЕЛЯ, ЖЕЛЕЗА И КОБАЛЬТА 

Обсудим теперь кратко результаты некоторых эксперименталь¬ 
ных и теоретических исследований ферромагнитных переходных 
металлов — №, Ре и Со. Общий обзор проблемы с различных точек 
зрения хорошо изложен в [27 —30]. 

Никель. Наиболее подробно исследован никель. Температура 
Кюри Т к (№) = 631 К [2]. Средний атомный магнитный момент 
для ферромагнитного состояния кристалла, выраженный в магне¬ 
тонах Бора [31]: р, (ІЧі) = 0,583 (данные для комнатной темпера¬ 
туры). Мотт одним из первых предположил, что наблюдаемые 
дробные магнитные моменты связаны с тем, что энергия Ферми 
лежит в пределах Зй-зоны. При этом энергия Ферми должна лежать 
на таких уровнях в N1, Ре и Со, где число дырок в Зй-зоне могло 
быть 0,6 в Ш, 1,7 в Со и 2,2 в Ре; это объясняет экспериментальные 
значения магнитного момента. Что может дать информация о сред¬ 
них атомных магнитных моментах для решения интересующего 
нас вопроса «локализованные — коллективизированные» электро¬ 
ны? Средний атомный магнитный момент является статической 
характеристикой. Теоретический расчет его обеспечивает хорошую 
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проверку для волновых функций электронов в металле [32, 33]. 
Степень локализации (или коллективизации) магнитоактивных 
электронов характеризуется лишь по тому, какой вклад в сред¬ 
ний момент вносят электроны, близкие ко дну или к верху зоны. 
Этого, однако, недостаточно, чтобы высказаться в пользу той 
или иной микроскопической модели (см. Введение). 

Вычисление температурной зависимости электронной тепло¬ 
емкости для ферромагнитных металлов впервые было выполнено 
Стонером [34] на основе зонной модели в приближении эффектив¬ 
ного поля и для параболической кривой плотности состояний. 
Этот расчет был усовершенствован Хантом [35], а впоследствии 
Шимицу и Терао [36] для случая никеля и железа. Низкотем¬ 
пературную магнитную теплоемкость никеля вычислил также 
Томпсон [37]. Было показано, что в области низких температур 
спин-волновая теплоемкость, вычисленная в приближении не¬ 
взаимодействующих спиновых волн, не вполне удовлетворительно 
согласуется с экспериментальными данными. Согласие удается 
улучшить [38] за счет подбора плотности состояний никеля, опи¬ 
рающегося на детальные теоретические расчеты зонной структу¬ 
ры [39—42]. Наилучшее согласие с экспериментом дают, однако, 
расчеты [36], выполненные на основе модели Стонера для никеля 
и железа. В этой модели существуют только стонеровские воз¬ 
буждения, которые в расчете Томпсона [37, 38] не принимаются 
во внимание. По мнению Шимицу и Терао [36] для правильного 
объяснения магнитной электронной теплоемкости при расчете 
должен учитываться вклад как от спин-волновых, так и от стоне- 
ровских возбуждений. 

За последние годы было выполнено довольно большое количе¬ 
ство теоретических расчетов зонной структуры ферромагнитного 
никеля [39—54]. На основе этих расчетов удалось удовлетворитель¬ 
но объяснить ферромагнитные энергии анизотропии никеля и же¬ 
леза [55], внутренние поля в никеле и железе и изомерный сдвиг 
железа [48], спин-волновые энергии железа, кобальта и никеля 
[52]; фактор спектроскопического расщепления и магнитомехани¬ 
ческое отношение никеля и железа [56] и целый ряд других харак¬ 
теристик. Наше понимание энергетической зонной структуры 
N1, Ре и Со было уточнено благодаря изучению рентгеновских 
спектров [57], измерению спиновой поляризации электронов, тун¬ 
нелирующих из ферромагнетиков [58—60], и фотоэмиссионным 
исследованиям [61—64]. Экспериментальное изучение фотоэмис- 
сионных электронов позволяет, в частности, определить зонное 
расщепление А. Для N1 наилучшая оценка по совокупности экс¬ 
периментальных и теоретических исследований А « 0,5 ±0,1 эВ. 
Фотоэмиссионные эксперименты, однако, не позволяют устано¬ 
вить вклад от различных типов возбуждений. Косвенные дока¬ 
зательства существования стонеровских возбуждений дают также 
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исследования эффекта де Гааза-ван Альфвена в переходных Зй- 
металлах [65, 66]. В работе [67] отмечалось, что точное описание 
комбинационного рассеяния света в переходных металлах требует 
учета добавочного рассеяния на стонеровских возбуждениях. 
При малых импульсах и нулевой температуре интенсивность рас¬ 
сеяния света на магнонах больше, чем на стонеровских возбужде¬ 
ниях. Однако вклад в интенсивность рассеяния света от магнонов 
и от стонеровских возбуждений по-разному зависит от специфики 
зонной структуры металла [68]. 

Железо. Температура Кюри 7\(Ге)=1043 К [2]. Средний 
атомный магнитный момент для ферромагнитного состояния 
р (Ре) = 2,177 [31]. Исследование температурной зависимости 
электронной теплоемкости железа [36] говорит о важности вклада 
от стонеровских возбуждений. Расчеты энергетической зонной 
структуры железа [69—73] позволили углубить наше понимание 
его магнитных свойств. Изучение поверхности Ферми железа 
[74—77] и никеля [74], спин-орбитальной связи и оптической 
проводимости железа [77] показывает очень хорошую примени¬ 
мость модели коллективизированных электронов. Исследования 
температурной зависимости обменного расщепления в железе [65] 
с помощью эффекта де Гааза — ван Альфвена дают указания на 
важность вклада от стонеровских возбуждений. 

Тем не менее единого мнения по поводу того, почему железо 
магнитно, все еще нет. Авторы работы [78] задают вопрос: должны 
ли магнитные свойства железа описываться на основе представле¬ 
ния о локализованных моментах или на основе зонного описания 
как в модели Стонера. В этой работе изучались сверхтонкие поля 
в сплавах ГеЗі и ГеМп. В целом информация о состоянии магнито¬ 
активных электронов в Ге, получаемая из измерения сверхтонких 
полей, проанализирована в обзоре [79]. Авторы работ [78, 79] 
убедительно показывают, что электроны в Ге отчасти сохраняют 
локализованный характер, и поэтому простая модель Стонера 
неприменима. Модель Хаббарда, учитывающая оба аспекта пове¬ 
дения электронов: коллективизированный и локализованный, 
значительно более адекватно описывает реальную ситуацию в ме¬ 
талле; в этой модели коллективизированный и локализованный 
аспекты поведения электронов носят дополнительный характер. 

Кобальт. Температура Кюри Т (Со) = 1403 К [2]. Средний 
атомный магнитный момент р (Со) = 1,707 [31]. Электронная 
энергетическая зонная структура кобальта была рассчитана в ра¬ 
ботах [80—83]. На основе этих расчетов дана простая интерпре¬ 
тация фотоэлектронных измерений в никеле и кобальте [84], а так¬ 
же интерпретация измерений сверхтонких полей [85]. Исследова¬ 
ния низкотемпературной электронной теплоемкости позволили 
определить коэффициент у [86]. Значение обменного расщепления 
кобальта больше, чем у никеля А (Со) == 0,91 эВ [80]. При этом 




382 а. Л. КУЗЕМСКИЙ 


показательно то, что отношение эффективного кулоновского взаи¬ 
модействия II к ширине й-зоны у всех трех металлов Ре, N1, Со 
приблизительно одинаково: И/го ж 0,14—0,16. 

В работе [87] зонная модель была использована для вычисле¬ 
ния температуры Кюри Со. Энергетическая зонная структура, 
волновые функции и плотность состояний находились из первых 
принципов. Полученное теоретическое значение Т к да 1370 ± 
± 200 К. Оценка для Со температуры перехода на основе модели 
локализованных спинов Гейзенберга дает завышенное значение, 
а именно Т к да 1870 К. Перейдем теперь к обсуждению данных 
о спектре магнитных возбуждений и намагниченности. 

5. СПЕКТР МАГНИТНЫХ ВОЗБУЖДЕНИИ Ш, Ре, Со 

Очень важная информация о состоянии магнитоактивных элек¬ 
тронов получена при исследовании спектра магнитных возбужде¬ 
ний. Наиболее интересные (с точки зрения сравнения локализо¬ 
ванной и зонной моделей) исследования были проведены с помощью 
рассеяния нейтронов. 

Рассеяние нейтронов в N1 детально исследовалось в ряде тео¬ 
ретических и экспериментальных работ [88—101]. Обзоры ранних 
нейтронных исследований спектра магнитных возбуждений в N1 
даны в [95, 100, 101]; более поздние результаты обсуждаются в об¬ 
зоре [30]. 

Метод рассеяния тепловых нейтронов позволяет непосредствен¬ 
но измерять закон рассеяния 5 (д, ю) в широком интервале значе¬ 
ний переданных энергий Йсо и импульсов д: 

^ = к 0 —к = (т/к)(ѵ 0 — ѵ); т/Н = 0,3956А -1 см _1 *мс; 

Иа = {П 2 /2т) {к 2 0 — к 2 ) = (1/2) т (ѵ 2 а — у 2 ); Я 2 /2ш = 2,072 МэВ-А" 2 . 

С помощью рассеяния тепловых нейтронов можно в принципе 
исследовать область значений д в интервале 0,1—10 А -1 и энер¬ 
гий Йсо в интервале 10— ІО 3 МэВ. 

При изучении рассеяния нейтронов в переходных металлах 
можно в первом приближении пренебречь вкладом орбитального 
момента [102]. Для N1 расчет показывает [96], что при больших 
передачах импульса и энергий значение орбитального вклада 
возрастает. Ферромагнитные металлы имеют хорошо выраженные 
пики неупругого рассеяния на спиновых волнах. Для N1, Ре, Со 
законТдисперсии 

йсо = І)д 2 |(1-Рд 2 + уд 4 ) (34) 

измерен с большой точностью. Результаты нейтронных измерений 
для № приведены на рис. 4. Заметим, что в ряде работ [103—106] 
спектр спин-волновых возбуждений, измеряемый с помощью рас- 
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сеяния нейтронов, интерпретировался на основе модели Гейзен¬ 
берга. Эксперименты проводились на N1 [104] и Ре [104—106]. 
Однако последующие исследования спектра спиновых волн в Ре 
[107—111] и Со [90, 100, 112, 115] показали, что коллективизирован¬ 
ная модель более соответствует действительности. Решающим аргу¬ 
ментом в пользу коллективизированной модели было косвенное 
указание на существование стонеровских возбуждений. 

Было обнаружено [99], что в N1 спин-волновая интенсивность 
рассеяния нейтронов вдоль направления [110] медленно спадает 



Рис. 5. Измеренная спин-волно¬ 
вая интенсивность в ІЧі при ком¬ 
натной температуре вдоль трех 
симметричных направлений [99] 


Рис. 6. Результаты теоретического 
расчета спектра возбуждений для 
N1 [98] в направлении [111] 



с увеличением энергии и при значении Е м 100 МэВ резко умень¬ 
шается более чем на порядок. Такое внезапное уменьшение интен¬ 
сивности было объяснено пересечением спин-волнового спектра 
и стонеровского континуума. Значение энергии, при котором про¬ 
исходило уменьшение интенсивности, оказалось различным для 
трех главных симметричных направлений (80 МэВ для направле¬ 
ния [111] и 110 МэВ для направления [100]) (рис. 5). 

Расчеты обобщенной восприимчивости Х +_ (д, со) и ее полюсов, 
проведенные на основе вычисления зонной энергетической струк¬ 
туры N1 [98] с учетом зависимости зонного спинового расщепления 
от квазиимпульса и многозонных эффектов, показали, что вдоль 
направления [111] область стонеровских возбуждений находится 
ниже, чем вдоль двух других направлений (рис. 6), что согласует¬ 
ся с экспериментом. Подобное же поведение спин-волновой интен¬ 
сивности наблюдается и в Ре [109] с тем отличием, что значения 
энергий спада Еж 90 МэВ для направления [100], Ем 95 МэВ 
для [110] иЕй 100 МэВ для [111]. 

Как для N1 [95, 99], так и для Ге [110] были проведены деталь¬ 
ные нейтронные измерения температурной зависимости спектра 
магнитных возбуждений. Ожидалось, что с изменением темпера- 
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туры стонеровский континуум будет сдвигаться, поскольку долж¬ 
но измениться значение зонного расщепления А, пропорциональ¬ 
ное намагниченности. Измерения в N1 [99] вдоль направления [111] 
показали, что точка пересечения спин-волнового и стонеровского 
Спектра практически не меняется в интервале температур от 4,2 
до 715 К (рис. 7). Более того спин-волновой спектр в этом темпе¬ 
ратурном интервале также изменяется незначительно. Коэффици¬ 
ент жесткости Б в законе дисперсии Йсо = Юср уменьшался от 



Рис. 8. Теоретический расчет за¬ 
висимости полного сечения рас¬ 
сеяния в N1 от квазиимпульса 



В 0,1 0,2 0,3 0,40 0,1 0,2 <}3 0,4 


555 МэВ А 2 при 4,2 К до 280 МэВА 2 при Т к = 631 К; при даль¬ 
нейшем увеличении температуры до 715 К коэффициент 7) прак¬ 
тически не изменялся. 

Эти результаты были получены с помощью трехосного спектро¬ 
метра. Техника малруглового рассеяния [107] дает другое значе¬ 
ние Б = 125 МэВ А 2 при Т = Т к . Однако независимо от этих 
расхождений выполняется условие Ъ ( Т ) -ѵ7) 0 п Ри Т -+Т к . 
Подобное же температурное поведение наблюдается и для Ре [110]. 
Для измерений использовались кристаллы 60 № и ( 54 Ре + 8і). 
Изотопы в0 № и 64 Ре имеют очень малую ядерную амплитуду рас¬ 
сеяния по сравнению с нормальными изотопами. При этом фонон¬ 
ное рассеяние здесь меньше и ядерное некогерентное рассеяние 
дает малый вклад. 

Таким образом, можно говорить о том, что эксперименты по 
неупругому рассеянию, нейтронов в переходных металлах дают 
важные косвенные указания на существование стонеровских воз¬ 
буждений. Подробные теоретические расчеты [98] поддерживают 
этот вывод. 
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Каковы перспективы непосредственного исследования стоне- 
ровских возбуждений? Важнейшим следствием расчетов, прове¬ 
денных в [98], является утверждение, что в области больших 
передач импульса спин-волновое сечение рассеяния всего в два- 
три раза больше стонеровского сечения рассеяния (рис. 8 ). Можно 
дать очень простую качественную интерпретацию этого обстоя¬ 
тельства 124, 25], которая основывается на результатах теорети¬ 
ческого расчета Томпсона [116]. Томпсон рассматривал модель 
Хаббарда, исходя из правила сумм (20): 

^ Іш Х + ~ (д, со) == (2я/Й) (/4 — п^) 


и использовал приближение хаотических фаз для расчета вос¬ 
приимчивости и приближение эффективной массы. Численным 
образом он показал, что имеет место равенство 

| йю Іпі Х^“ (д — 0,9д макс , = 


— | с?со Іш Хі„. (д = 0,9д макс , со) = 

— 2я А (л _1\ __2я_ А 

— Н 2 )~ К 211 ’ 


(35) 


При этом Томпсоном было получено, что для сильного ферромагне¬ 
тика (Д > г } , что соответствует N 1 ) свободная восприимчивость 


Іт Х^" (д, со) = 


{ 


0 , 

_ 3 _ 

4 




\У\> 1; 

Ш<і, 


гд еу = (А' — со)/бі57. А' = А Н\ 2 /2т* и б Е = Н\Ц1т*’, т* — 
эффективная масса электрона. Таким образом, ширина диффузно¬ 
го пика, определяемая Іш Х„ _ (д, со), будет 8Е = (д/Ау) 8 у. 

Оценим величину б Е. Следуя Маттису [117], примем, что 
д ма кс < 0,75Тогда для ^ = 0,9д макс найдем, что б Е «: 0,67бу. 
Точные расчеты 8 у для металлического N1 [50] показывают, что 

= 7,6 эВ. Однако для оценки Іш Х + 0 _ в рамках расчета Томп¬ 
сона, чтобы не нарушать «самосогласованности» приближений, 
воспользуемся оценками 8 у для так называемой модели эквива¬ 
лентного свободного электронного газа в парамагнитной области 
[95]. В рамках этой модели, отвечающей приближениям работы 
[116], 8 у = 0,4 эВ. Таким образом, ширина стонеровского диф¬ 
фузного пика б Е — 0,67бу = 0,67-0,4 = 0,27 эВ. Ширину спин- 
волнового пика при д = 0,9д макс можно принять равной 0,1 эВ 
[94, 95, 99, 100]. Поскольку согласно формуле (35) площади обоих 


і — 01068 
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пиков равны для д = 0,9д мако , следовательно, амплитуда стоне- 
ровского пика будет примерно в 3 раза меньше амплитуды спин- 
волнового пика. Таким образом, проведенные качественные оценки 
находятся в согласии с точными численными расчетами (рис. 8) 
работы [98]. 

Измерения спектра магнитных возбуждений в N1 [99] и Ре 
[109, 110] были выполнены на стационарных реакторах при 
помощи трехосного кристаллического спектрометра. Зонное рас¬ 
щепление А N1, Ре и Со имеет наименьшее значение для N1 А да 
« 0,5 эВ. Поэтому для измерения стонеровских возбуждений 
в N1 необходим поток нейтронов с энергией в интервале 0,1 — 
1 эВ. Напомним, что точке окончания спин-волнового спектра 
вдоль направления [111] соответствует энергия 80 МэВ, поэтому 
интересующая нас область передачи энергий заключена в интерва¬ 
ле 0,08—0,5 эВ. 

Ясно, что обычная техника мало подходит для измерения воз¬ 
буждений столь высокой энергии даже при наличии горячего 
источника. Многообещающим методом является метод времени 
пролета. Импульсные источники нейтронов имеют необходимую 
мощность для обеспечения достаточной светосилы, а малая шири¬ 
на импульса позволяет работать при хорошем разрешении и низ¬ 
ком уровне фона. Согласно оценкам работы [118] при Е 0 = 0,5 эВ 
и А Е 0 — 0,0045 эВ, времени вспышки т = 3 мкс и размерах замед¬ 
лителя 20 X 20 см поток нейтронов, ожидающийся на расстоянии 
30 м от ИБР-2, можно считать равным 

/ = 2,5-10 4 нейтр./(см-с 2 ). 

В работе [24] были проведены оценки скорости счета установ¬ 
ки при рассеянии такого потока медленных нейтронов на стоне¬ 
ровских возбуждениях. При этом использовались те оценки для 
стонеровского сечения рассеяния, которые обсуждались выше. 
Был получен следующий результат: ожидаемое число отсчетов 
за время измерения I = 24 ч N та 200—250 нейтр./сут. Аналогич¬ 
ные оценки для мощного импульсного источника нейтронов 2ШС 
[119] дают N « 160 нейтр./сут. Таким образом, можно говорить, 
что порядок величины N « ІО 2 нейтр./сут. 

Приведенные оценки сверху соответствуют высокой разрешаю¬ 
щей способности, позволяющей детально исследовать свойства 
стонеровских возбуждений. Если же целью эксперимента является 
простое обнаружение их существования, то достаточно работать 
с разрешением в десять раз хуже. При этом интенсивность вырастет 
на порядок и скорость счета повысится. 

На имеющихся в настоящее время установках непосредственно 
измерять стонеровские возбуждения в №, Ее, Со довольно трудно 
из-за большого значения зонного расщепления. Поэтому экспе¬ 
риментаторы стараются найти зонные ферромагнетики с неболь- 
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шим значением Д. В работе [120] для исследования магнитных 
возбуждений выбрано интерметаллическое соединение МпЗі, кото¬ 
рое является слабым зонным ферромагнетиком. Ферромагнитное 
состояние в МиЗі индуцируется внешним магнитным полем с нап¬ 
ряженностью выше 6 кЭ. Эксперимент по рассеянию нейтронов 
был проведен с помощью трехосного кристаллического спектро¬ 
метра. 

В магнитном поле напряженностью 10 кЭ при температуре 5 К 
наблюдались хорошо определенные спйн-волновые возбуждения 
ниже 2,5 МэВ. Закон дисперсии был найден в виде Йсо 9 (МэВ) = 
= 0,13+ 52# 2 (А 2 ) для направления [100]. Выше 3 МэВ наблю¬ 
далось значительное уширение спин-волнового пика, что связы¬ 
валось авторами с пересечением стонеровского континуума. Таким 
образом, коэффициент жесткости в МпЗі на порядок меньше, чем 
в никеле. Зонное расщепление определить не удалось; была найде¬ 
на только энергия пересечения стонеровского континуума. Для 
направления [100] эта энергия равна 2,6 МэВ. 

Рассеяние нейтронов в Ре 3 А1 [121] также дает косвенное под¬ 
тверждение существования стонеровских возбуждений. Для нап¬ 
равления [111] энергия пересечения стонеровского континуума 
около 12 МэВ. В интерметаллическом соединении Р6 3 Ре [122] 
кривая дисперсии магнонов в направлении [001] пересекает стоне- 
ровский континуум при Еж 40 МэВ. 

Итак, исследования спектра магнитных возбуждений ферро¬ 
магнитных переходных металлов и некоторых их соединений дают 
очень веские аргументы в пользу того, что стонеровские возбуж¬ 
дения существуют и с развитием техники эксперимента их мож¬ 
но будет наблюдать непосредственно. 

6. ВКЛАД СПИНОВЫХ ВОЛН И СТОНЕРОВСКИХ ВОЗБУЖДЕНИЙ 
В НАМАГНИЧЕННОСТЬ 

Большой интерес, с точки зрения сравнения локализованной 
и коллективизированной моделей, представляет изучение темпера¬ 
турной зависимости намагниченности. К сожалению, до сих пор 
нет достаточной ясности в этом вопросе. 

Вольфарт [27] проанализировал магнитные, тепловые и магни¬ 
тоупругие свойства сильных и слабых зонных ферромагнетиков, 
исходя из простых термодинамических соображений. По его мне¬ 
нию, свойства сильных зонных ферромагнетиков, в частности N1, 
ниже точки Кюри в основном хорошо описываются на основе 
простой спин-волновой теории. Стонеровские возбуждения имеют 
существенное значение в описании слабых ферромагнетиков. 

Более детальный анализ Шимицу [28] и в особенности Реди [123] 
показывает, что ситуация гораздо более сложная, чем это следует 
из работы [27]. 
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Реди [123] анализирует совокупность экспериментальных дан¬ 
ных по рассеянию нейтронов,- намагниченности, ЯМР и пытается 
сделать заключение об относительном вкладе спиновых волн 
и стонеровских возбуждений в температурное Поведение намагни¬ 
ченности М (Т). 

В самом деле, температурная зависимость намагниченности 
существенно зависит от модели, которую будем использовать для 
описания. Для простого кубического ферромагнетика Гейзенберга 
температурная зависимость намагниченности в спин-волновом 
приближении имеет вид 

АМ а .^М(0)-М(Т)~Г(3/2Т)(е\ів/М(0)) ( 1 ^ т?г ) 3/2 _ 

- О (Т 5,2 ) = аТ 3/2 (1+ЬТ +...). (36) 

Впервые коэффициенты а при члене Т 3/і для N1 и Ре были 
найдены в работе [124]. В зонном ферромагнетике, помимо спино¬ 
вых волн, существуют также стонеровские возбуждения, которые 
также дают вклад в термодинамическое поведение. Поэтому в об¬ 
щем случае температурная зависимость ДМ ( Т ) должна опреде¬ 
ляться вкладом от спиновых волн и стонеровских возбуждений, 
а также от их взаимодействия. 

В некоторой области значений параметров системы ДМ ( Т ) 
можно приближенно представить в виде суперпозиции этих двух 
вкладов: 

ДМ (Т) ~ Д М в . ѵ . + ДМа.р.. (37) 

Температурная зависимость ДМ 3 . Р .(Г) отличается для слабого 
и сильного ферромагнетиков (см. [123]): 

ДМ 8 .р. (Т) = 

{ п 00 [21 ( Т)/п] ехр (— Д/М 1 ) — сильный ферромагнетик; 
о 00 АТ 2 —слабый ферромагнетик. 

Более точный расчет Херринга [125] с учетом взаимодействия 
спиновых волн с одночастичными возбуждениями дает для слабого 
ферромагнетика следующее выражение: 

АМ &Л> .{Т) = АП-СТ Ь ' 2 . (39) 

Таким образом, исследование температурной зависимости на¬ 
магниченности в принципе может дать важные аргументы в пользу 
справедливости той или иной модели. Дело, однако, осложняется 
тем, что для надежных заключений необходима очень большая 
точность измерений, поскольку ДМ ( Т) очень мала в низкотем¬ 
пературной области. 

В последние годы были проведены подробные теоретические 
и экспериментальные исследования температурной зависимости 
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намагниченности N1 [124, 126—130] и Ре [124, 128, 131—134]. 
Анализ данных по рассеянию нейтронов, измерению намагничен¬ 
ности и ЯМР показывает, что N1 и Ре ведут себя по-разному. 
Низкотемпературная намагниченность Ре хорошо описывается 
формулой 

Ш {Т) = В 0 Т 3/2 , (40) 

где В 0 — константа (в согласии со спин-волновой теорией). 

Как следует из (36), измеряя с помощью неупругого рассеяния 
нейтронов спин-волновую жесткость 7), можно найти коэффициент 
В 0 . Для Ре все три значения В 0 , получаемые с помощью неупруго¬ 
го рассеяния нейтронов, измерения намагниченности, ЯМР, хоро¬ 
шо согласуются между собой. 

При более высоких температурах коэффициент В в (40) зависит 
от температуры. Для модели Хаббарда расчет коэффициента спин- 
волновой жесткости в приближении Хартри — Фока приводит 
к следующей температурной зависимости [4]: 

О(Т) = О 0 — О і Т 2 -О 2 Т 5 ' 2 . (41) 

Вклад в И (Т), пропорциональный Т 2 , возникает за счет учета 
электрон-магнонного взаимодействия. Стрингфеллоу [107] с помо¬ 
щью техники малоуглового рассеяния измерил 77 ( Т) для Ш и Ре. 
Было найдено, что 7^ = 0 для Ш. Для Ре 7Д Ф 0 в низкотемпера¬ 
турной области и зависимость 7) от Т 2 существенна. Таким образом, 
температурная зависимость спин-волновой жесткости для Ш 
имеет вид 

7? = 7?о-7) 1 Г 5/2 . (42) 

Согласно Альфреду наилучшие значения 7) для Ш: 7) (Т = 4,2 К) =» 
= 555 МэВ- А 2 и 7) (Т = 295 К) = 455 МэВ-А 2 . На основе этих 
значений, полученных с помощью рассеяния нейтронов, можно 
определить коэффициент В температурной зависимости намагни¬ 
ченности (40). 

Оказалось, что для Ш величина В , полученная из измерений 
намагниченности и ЯМР, примерно на 30—40% меньше, чем то 
значение, которое дают нейтронные эксперименты. Альдред [129] 
высказал предположение, что такое резкое различие в поведении 
N1 и Ре можно объяснить тем, что низкотемпературное поведение 
Ре удовлетворительно описывается спин-волновой картиной, но 
для № необходимо учитывать вклад от стонеровских возбуждений. 
Поэтому Альдред рассмотрел добавочный стонеровский вклад 
в намагниченность Ш: 

ДЛ7 8 .р. (7 1 ) = АМ (Т) — АМ (ГД.-йг. • (43) 

Измеренную температурную зависимость АМ ( Т) он попытался 
описать на основе выражения (38). К сожалению, точность изме- 
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рений [129] не позволила отдать предпочтение какой-либо из воз¬ 
можностей в (38). Обе зависимости ДМ 8 . Р . не вполне удовлетвори¬ 
тельно ложатся на экспериментальную кривую, хотя зависимость 
АТ 2 все же несколько лучше приближает ее. 

Реди [130] определил вклад ДМ 8 . р . для N1 с помощью техники 
ЯМР. Его результат имеет вид 

дм 8 .р. = 1,68 • 10" в Г 3/2 + 3,22• 10- 7 Г 2 . (44) 

Реди показал, что член, пропорциональный Т 2 , не зависит от Б ( Т) 
во всей области экспериментальных значений спин-волновой жест¬ 
кости. Поэтому, считает Реди [123, 130], можно связать зависи¬ 
мость Т 2 в (44) с вкладом стонеровских возбуждений. Анализ 
Реди показывает, что выражение (39) хуже согласуется с экспери¬ 
ментом, чем простая зависимость Т а . 

Несомненно, результаты работ [123—134] имеют очень большое 
значение, хотя не вносят полной ясности. Из (44) следует, что 
величина Д /к Ф 0 очень мала, поскольку Ш является сильным 
ферромагнетиком [135]. Измерения намагниченности дают А/к = 
= 162 К [130]. Реди полагает, что это значение должно быть, по 
крайней мере, на порядок меньше, чтобы объяснить наблюдаемое 
расхождение с нейтронными данными. 

Таким образом, для окончательного суждения о вкладе стоне¬ 
ровских и спин-волновых возбуждений нужны дополнительные 
очень точные эксперименты. Ясно, однако, что для N1 только 
спин-волнового описания недостаточно. Надеяться, что будет 
достигнута необходимая точность измерений температурной зави¬ 
симости намагниченности, достаточная для однозначных выводов, 
трудно, поскольку оба вклада входят аддитивным образом. Вот 
почему мы считаем, что только метод рассеяния медленных нейтро¬ 
нов может дать надежную и прямую информацию о стонеровских 
возбуждениях [136]. 

7. НЕКОТОРЫЕ ОБОБЩЕНИЯ МОДЕЛИ ХАББАРДА 

При обсуждении экспериментальных данных по рассеянию 
медленных нейтронов в ферромагнетиках необходимо учитывать 
приближенный характер теоретических предсказаний, обсуждав¬ 
шихся выше. До сих пор речь шла о спектре возбуждений однозон¬ 
ной модели Хаббарда. На самом деле в реальном металле ситуа¬ 
ция может быть сложнее. Уже для двухзонной модели картина 
неупругого рассеяния нейтронов качественно и количественно мо¬ 
дифицируется; при этом стонеровское сечение рассеяния может 
стать заметным и для не слишком больших д < д ма _ КС . 

Спектр магнитных возбуждений двухзонной модели Хаббарда 
[137] показан на рис. 9. В отличие от однозонной модели возника¬ 
ют четыре квазистонеровских континуума, связанные с перехода- 
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ми внутри одной зоны и между зонами, а также ветвь оптических 
спиновых волн и так называемая межзонная спин-волновая ветвь. 
Таким образом, модель коллективизированных электронов облада¬ 
ет весьма богатым спектром возбуждений в «оптической» области. 

1. Важным обобщением однозонной модели Хаббарда является 
модель Хаббарда с (в — гі)-гибридизацией, описывающей прямое 
рассеяние электронов к- и й-зон [138]: 

Я= 2 Міо^о + 4 2 «іо«і-а + 

+ 2 Е (&) а Асг а йсг + 2 (^ 2 Ь^Я» а Ьо + Э.С.). ( 45 ) 

к іка 

Существование широкой зоны 5-электронов неявно учитывается 
и в обычной однозонной модели Хаббарда при реалистических 
оценках параметра кулоновского оттал¬ 
кивания. Поскольку скорость пере¬ 
скока й-электронов с атома на атом зна¬ 
чительно меньше характерной скорости 
электронов проводимости, последние мо¬ 
гут эффективно коррелировать с й-эле- 
ктронами и экранировать их поля. Для 
переходных Зй-металлов исследование 
энергетической зонной структуры пока¬ 
зывает важную роль процессов 5 — д,- 
гибридизации. Займан показал [139], 
что для благородных и переходных 
металлов в рамках применимости метода 
ККР й-зона фактически является резо¬ 
нансом в широкой зр-зоне. Это позво¬ 
ляет предположить, что уширение атом¬ 
ных ^-уровней может проистекать из-за 
5 — й-гибридизации, конкурируя с пря¬ 
мым перекрытием волновых функций 
^-электронов. Эта точка зрения полу¬ 
чила подтверждение в работах [140—143]. Впоследствии прямые 
расчеты на основе интерполяционных схем [144] показали, что 
интегралы перекрытия волновых функций й-электронов могут быть 
меньше, чем интегралы перекрытия волновых функций 5 - и й-эле- 
ктронов, хотя отличие не слишком значительно. 

Таким образом, в ряду других обобщений модели Хаббарда, 
более соответствующих реальной ситуации в чистых переходных 
металлах, модель с 5 — ^-гибридизацией занимает одно из глав¬ 
ных мест. Одночастичные свойства различных вариантов этой 
модели были исследованы в [145—148]. Спин-волновой полюс 
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Рис. 9.] Спектр возбужде¬ 
ний двухзонной модели 
Хаббарда, рассчитанный в 
[137] 
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обобщенной спиновой восприимчивости для модели (45) вычислял¬ 
ся в [149] в атомном пределе для й-электронов и для очень сильной 
5 — й-гибридизации. Однако эти приближения совершенно непри¬ 
менимы для реальных переходных металлов. Подробная теория 
рассеяния медленных нейтронов для модели Хаббарда с $ — й- 
гибридизацией приведена в [150]; данный расчет является прямым 
обобщением работы Изюяма, Кима и Кубо [26]. Они показали, что 
для однозонной модели Хаббарда, пренебрегая перекрытием вол¬ 
новых функций на разных узлах, можно выразить обобщенную 
спиновую поперечную восприимчивость Х + ~ (с{, со) системы через 
хартри-фоковскую восприимчивость Х+" {д , со) в виде (24). 

Хартри-фоковская восприимчивость системы (45) вычисляется 
с помощью простого (и, ^-преобразования линеаризованного 
гамильтониана 

Н НР = 2 Ё а (к) Йй а Й йсг 4~ 2 ^к а ко^ко + 2 (Ук а ка^каЛ- Э.С.) (46) 

Аа ко ко 

к следующему виду: 

дгНР_ ^ {со 1йо ад ст а йо + ЮгйоРАстРйа}- (47) 

ко 

Здесь 

< 02 *а (к) + е*) ± Ѵ(Ё а (к)-е к )* + 4\Ѵ к \>; (48) 

І л да=г(і)+І2ѵ 


Выражение для мнимой части хартри-фоковской восприимчивости 
записывается в форме [150]: 


ІщХ НІ? (д, (й)== — л(^р в ) 2 -^ ^ 'п“ +д |) X 

А 

X 6 (ЙСО + ООій+д^ — С02А^) + 

+ (»А| — »А+в|) & (ЙСО + (02А+д| — ®2А|) + 

+ и\+і\ ѵ 1\ ( К А| — И А+^) б (ЙСО + С01А+д^ — (І>2к\) + 

+ ѵ к+9 \иІ^(п к ^ — п к+а ^)& (йсо + соіа+,4 — ©2а|)}- (49) 


Таким образом, стонеровское сечение рассеяния медленных ней¬ 
тронов 


(_&о_\ _ / уе 2 

V йййЕ /а.р. V т е с 2 

Х{(я (со) -Ь1) ІтХ нг 


■) Чі'мі’т-^+йтіет-х 

( — д, со) + п (— со) Іт Х НР (д, со)} 


(50) 
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будет описывать рассеяние на четырех квазистонеровских конти¬ 
нуумах. 

Чтобы рассчитать обобщенную спиновую восприимчивость от 
поперечных спиновых компонент для модели (45), необходимо 
вычислять двухвременную функцию Грина «Ѳ й (?) = «+,^ м |Я» в . 
Цепочка уравнений для функций Грина сводится к замкнутой 
системе четырех уравнений с помощью следующего приближения 
хаотических фаз: 

[Ѳ* (?), Н а ] » (Е (к) -Е(к + ?)) Ѳ к (?) + 

+ АѲ Й (?) — 2 ( п *+<г1 - п м) Ѳ Р (5)'» 

V 

Ык+^\ а к\1 Нд] « Е (А-)- ?) (І к +д\ а к\ 2 п Ѵ\^к+а\ а к\ + 

Р 

+ 7Г {(Ік\ а к\) 2 ^Р + 9І^р1" 

В результате выражение для поперечной спиновой восприим¬ 
чивости Х -+ (?, со) запишется в виде, аналогичном (24): 

Х- + (?, со) = X й ' р (?, со) {1 - 1Л{ ё Ы 2 X я * (?, со)}" 1 . (51) 

Здесь восприимчивость в приближении Хартри — Фока Х нр (?, со) 
имеет сложный вид (см. [150]) и в пределе V Ік -ѵ0 переходит 
в обычное выражение (25). Поскольку гамильтониан (45) остается 
ротационно-инвариантным, то в силу общей теории Эдвардса (17) 
восприимчивость (51) должна иметь спин-волновой полюс. Наличие 
спин-волнового полюса у восприимчивости Х -+ (?, со) показано 
в [150]. При этом имеем 

ІтХ- + ( ? -*»0, со-»-0)-я(^ц в ) 2 Ь/Е7б(Йсо — ЙЙ,), (52) 

где Ъ — некоторая константа, определяемая численным образом. 
Сечение неупругого рассеяния нейтронов на спиновых волнах 
записывается в данном случае совершенно аналогично (32). 

2. В зонной модели магнетизма, как и в локализованной моде¬ 
ли, большой интерес представляет учет в явном виде электрон- 
фононного взаимодействия. Благодаря этому взаимодействию 
происходит неупругая передача энергии от системы магнитных 
электронов колебаниям решетки; при детальном описании рас¬ 
сеяния нейтронов это обстоятельство имеет существенное значе¬ 
ние [151]. 

Электрон-фононное взаимодействие определяется изменением 
потенциалов ионов решетки ѵ (К) при учете их тепловых колеба¬ 
ний: 

Н е - Р п = 2 {ѵ (В, - г г ) - <н (Кг - г ,)>,*}, 


( 53 ) 
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где координаты ионов К* = 1 + иг, 1 — равновесные положения; 
и; — тепловые смещения; г г — координаты электронов. В [151] 
блоховские функции определены по усредненному по колебаниям 
решетки потенциалу ионов и (1 — г) = (и (К — г) ) рН 

+ 0}фк(г) = ЕкФк(г). (54) 


Это определение отличается от обычного [см. (1) — (3)] и позволяет 
учитывать упругие процессы электрон-фононного рассеяния. 
В представлении вторичного квантования гамильтониан однозон¬ 
ной модели Хаббарда с учетом электрон-фононного взаимодействия 
записывается в виде [151]: 


н — Не -[- Я рЛ 4- Не-рП, 

(55) 


(55а) 

Н е -рН — 2 ^Р. р+Я а РО а р+ча X 


X 2 ех Р (ІЧ 1 ) ( ех Р (І4 ц г) — ( ех Р (ІЧЧ/))) • 

(556) 


Здесь Н е — однозонный гамильтониан Хаббарда (И); ѵ = ІЛЫ — 
энергия кулоновской корреляции в одном узле; V = — 

объем системы; Ь, и Ъ д — операторы рождения и уничтожения 
фононов с квазиимпульсом ц, поляризацией у и энергией а> д] ; 
4 = (д, у). Матричный элемент электрон-ионного взаимодействия 
имеет вид 

Ѵ р.ѵ+ч = ѵ(я) 5 #пр*(г)ехр( — іцг) Ф Р+Ч (г) = ѵ (?) Р (?, р + д ) » Ѵ д , 

где ѵ (<у) = 1/Ѵ 0 | сРг ехр (— іцг) ѵ (г) — фурье-компонента элек¬ 
трон-ионного потенциала. Потенциал ѵ (К; — г) имеет смысл 
эффективного (экранированного электронами проводимости) по¬ 
тенциала взаимодействия магнитных электронов и ионов. 

Н. М. Планида и Л. С. Смирнов [151] для системы с гамильто¬ 
нианом (55) вычислили сечение магнитовибрационного рассеяния 
при учете однофононного перехода: 


(и-).'”- <*> Т-Н І- 2 "' <»» тлг 2 * 

т РІ 

< [га(«> р у) +-|-±-і-]8(ц±р— 2ят)б(ш ± 10 ^). (56) 


Здесь 


РтЛ*) = (Ѵг 0 )М1-еі)хІ^( & )І , І^І а т[ ] 2 


(57) 
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— эффективное сечение магнитовибрационного рассеяния. Как 
видно, сечение рассеяния (56) имеет вид обычного однофононного 
ядерного рассеяния. Однако эффективное сечение (57) помимо 
•обычных для упругого магнитного рассеяния множителей 
((1 — еі) | Р (к) | 2 ) и квадрата намагниченности (ге* — гс*) 2 ~ 
^-(я^ —Я|) 2 , как и в модели локализованных спинов, допол¬ 
нительно содержит квадрат матричного элемента электрон-фонон- 
ного взаимодействия | V ц | 2 , что позволяет оценить его по магнито- 
вибрационному рассеянию нейтронов в металлах. 

Электрон-фононное взаимодействие в модели Хаббарда рас¬ 
сматривалось также в [152—156]. Морковский [153] вычислил вре¬ 
мя релаксации магнонов вследствие взаимодействия с фононами 
во втором порядке теории возмущений. Джордж [154] с помощью 
метода двувременных температурных функций Грина получил 
выражение для ренормированной спин-волновой жесткости в при¬ 
ближении хаотических фаз. В работе Ямады [156] вычислялись 
время релаксации магнонов и ренормировка спин-волновой жест¬ 
кости с помощью мацубаровских функций Грина в приближении 
хаотических фаз. Существенно, что в отличие от других авторов 
[154, 155] расчет Ямады [156] дает выражение для б Ь ~ Г 4 , кото¬ 
рое согласуется с феноменологическим результатом Изюяма и Кубо 
[152]. Ямада [156] также дал численные оценки указанных величин 
для N1 и Ре. 

3. С точки зрения рассеяния медленных нейтронов большой 
интерес представляет спектр магнитных возбуждений так назы¬ 
ваемой модели Зенера [157]. Гамильтониан модели Зенера запи¬ 
сывается в виде 

Я = 2 і і} а\ 0 а ]а + [7/2 2 п іа п і _ а -2^ 2 о г 8 г , (58) 

На іа і 

где а+ — а^ац] о~ = ст| = 1/2 (аДа^ — ацац) 

— спиновые операторы коллективизированных электронов; 8 г — 
оператор локализованного спина в узле і, отвечающий правилу 
Хунда. Перепишем гамильтониан (58) в ^-представлении 

Н = 2 г к а каО'ка + 77 ІѴ -1 2 ^к+я\ а к^ а к'~я\ а к’\ — 

— ^N 1/2 2 1^к+я^ а к1^я 0,к+я\ а к\^я~^(. а к+я\ а к \ — Як+яІ а ы) *$«] • (59) 

кя 

Обобщенная поперечная спиновая восприимчивость системы 
с гамильтонианом (59) 

Х + ~(7» ш )— | ехр (Ш) йі 


( 60 ) 
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вычисляется в приближении хаотических фаз, как и для обыч¬ 
ной модели Хаббарда (11). Полученный результат имеет 
вид [157] 

Х +- (д, ю) = Хо(д, ю)/{1 — \У— 2/ 2 <5 г >/(ш — .Тпх)] х 0 (д, ю)}, (61) 
где 

(62) 

— восприимчивость коллективизированных электронов в прибли¬ 
жении Хартри — Фока; Д = СТ/гж + 2/ (5 2 ); пх — А -1 X 
X 2( п м — п йі ). Аналогично находится поперечная восприим¬ 
чивость для локализованных спинов: 

<*5 + (д) 15" (— д)» ш = 2 <*5 2 >/{со — /пж + 

+ 2Л (5 2 > Хо (д, (о) [1 - Пхо (д, т)]“*>. (63) 

Для д, со->-0 найдем спектр спин-волновых возбуждений систе¬ 
мы, который определяется полюсами восприимчивости (61) и (63): 
й(о ас = Пд 2 —акустическая ветвь; ) 

Н(й 0р = Е 0р +И (ѴЕ 0р //А — 1) д 2 — оптическая ветвь. } 

Здесь Е ор = Хпх 2/ (5 2 ). Физически этот результат соответ¬ 
ствует тому, что в системе имеются неэквивалентные спиновые 
подсистемы коллективизированных и локализованных спинов. 
Помимо этих двух типов возбуждений имеется еще континуум 
стонеровских возбуждений 

/і(і^ = (би — Б]{-(^) -[- А а (65) 

который определяется полюсами восприимчивости (62). Схема 
спектра магнитных возбуждений модели Зенера дана на рис. 10. 
Таким образом, в отличие от однозонной модели Хаббарда (11) 
спектр модели Зенера содержит дополнительную ветвь оптических 
спиновых волн. Наличие этой оптической спиновой волны связано 
с тем, что наша система состоит из двух подсистем, как и для двух¬ 
зонной модели Хаббарда [137, 158]. Для того чтобы более ясно 
увидеть это, положим в гамильтониане (59) V = 0. Тогда придем 
к гамильтониану (я — д)- модели: 

Н = 2 Ча-коЧо — 2 ех Р (— І4 к г) X 

ко кяг 

X [<4+910/4*57 + Ок+я\а к ^8+ + (а*+д| а м- — <4+д|Ф4) *5|]. (66) 

Гамильтониан (66) широко используется для описания (я — і) 
обменного взаимодействия между спином примеси и спинами элек¬ 
тронов проводимости в немагнитных металлах; (см. [159] и цитиро¬ 
ванную там литературу). 
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Рассмотрим вычисление функций Грина ((5 _ (ц) | 5 + (— д)» ш . 
Для мнимой части восприимчивости в приближении хаотических 
фаз получим [25] 

— Іт %~ + (ц, со) = 256 (со- сОдс (ц)) + пхб (со — <а 0р (ч)) + 

+ (25 + пх) | й 3 ?[6(со — е к+аі — е й |) — 6 (со — е** + е й _ д! )]. 

(67) 


В приближении эффективной массы е к = к 2 /2т* при д - 
найдем, что 

п [1 — (2/3) 8/ (х/І8)] 


с (?) = ^і? 2 ! = 


®о Р (?) = / (25 + пх) + Б г д^\ Э 2 = 


2т* (28 + пх) ’ 

8 [1 + (4/3) е,/(/п)] 


г*х (28+пх) 


( 68 ) 


Таким образом, модель Хаббарда допускает большое много¬ 
образие обобщений, что позволяет учитывать целый ряд физически 
важных взаимодействий в реальном пере¬ 
ходном металле. Это особенно существенно 
при сравнении экспериментальных данных 
по неупругому рассеянию нейтронов в 
металле с результатами теоретических 
расчетов. Как следует из нашего обсуж¬ 
дения, детальное исследование спектра 
магнитных возбуждений в оптическом 
и акустическом диапазонах представляет 
очень большой интерес с точки зрения 
сравнения локализованной и зонной моде¬ 
лей магнетизма. Важно еще раз подчерк¬ 
нуть, что эти два подхода, локализован¬ 
ный и зонный, не противоречат друг 
Другу. Они, скорее, являются «дополни¬ 
тельными» аспектами квантовостатического 
описания состояния магнитоактивных эле¬ 
ктронов в реальном металле. Акустичес¬ 
кая часть спектра магнитных возбужде¬ 
ний является отражением определенной 
локализованности спинов, а оптическая — 

их делокализации. Эту глубинную «дополнительность» хорошо 
отражает модель Зенера (58), являющаяся одной из наиболее 
интересных моделей магнетизма. Разнообразие моделей отражает 
разнообразие интересующих нас [аспектов магнитного поведения 
реальных магнетиков. 



Рис. 10. Схематический 
вид спектра возбужде¬ 
ний модифицированной 
модели Зенера [157] 
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8. СПЛАВЫ МАГНИТНЫХ ПЕРЕХОДНЫХ МЕТАЛЛОВ 

В последние годы интенсивно изучали электронную структуру 
и разнообразные физические свойства сплавов переходных метал¬ 
лов. Для изучения магнитных свойств сплавов переходных метал¬ 
лов очень полезным оказался метод рассеяния медленных нейтро¬ 
нов. Исследование упругого и неупругого рассеяния медленных 
нейтронов в сплавах позволяет получить уникальную информацию 
о магнитных моментах и форм-факторах, а также об изменении 
спин-волновой жесткости. 

Необходимо отметить, что нейтронные исследования распре¬ 
деления магнитного момента в магнитных сплавах и изменения 
спин-волновой жесткости во многом стимулировали развитие 
современных методов расчета электронной структуры неупорядо¬ 
ченных сплавов, которые чрезвычайно полезны для решения мно¬ 
гих задач физики твердого тела. К ним относится широко теперь 
известный метод когерентного потенциала [160]. 

Модель Хаббарда оказалась очень полезной для описания мно¬ 
гих электронных и магнитных свойств сплавов переходных метал¬ 
лов и успешно применяется в большом числе работ. При описании 
неупорядоченных сплавов с помощью модели Хаббарда вводятся 
случайные параметры, поэтому говорят о модели Хаббарда со 
случайными параметрами. 

Перейдем к ее описанию. Предполагается, что взаимодействие 
электронов в бинарном неупорядоченном сплаве из двух магнит¬ 
ных компонент А С В ]_ с описывается следующим модельным гамиль¬ 
тонианом: 

Н = 2 Іі]йіаО-]а 4" 2 е і га іа4- 2 ^ і п іа п і-<у (69) 


Здесь, как и в (И), а іа , а* а — операторы уничтожения и рождения 
электронов Ванье в узле і со спином а. Считается, что интегралы 
перескока Іц одинаковы для обоих сортов атомов А и В, т. е. 
1 А А. = % в = 1 1} \ зонная структура чистых компонент А и В 
в отсутствие кулоновского взаимодействия одинаковая. Величины 
е,- и Ѵі — одночастичный потенциал и внутриатомное кулоновское 
взаимодействие соответственно: 



Для неупорядоченного сплава величины г г и 17 ( принимают слу¬ 
чайные значения в зависимости от того, заполнен ли узел атомом 
А или В. 

Гамильтониан (69) исследовали многие авторы в различных 
предельных случаях. Если предположим, что какая-либо из ком- 
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понент сплава (например, В) состоит из немагнитных атомов, то 
можно положить параметр Ѵ в = 0. Этот случай соответствует 
модели Вольффа [161, 162]. Если положим = е® в (69), полу¬ 
чим модельный гамильтониан, который рядом авторов [163, 164] 
был использован для теоретического описания сплава Рй — N1. 
Случай, когда — С/?, рассмотрен Лютером и Фульде [165] 
для анализа рассеяния парамагнонов на примесях; Ямада и Шими- 
цу [166] рассчитали спин-волновой спектр. Мория [167] детально 
исследовал электронную структуру вблизи магнитной примеси 
( V і' Ф 0) в немагнитной матрице (V® = 0) и рассчитал целый 
ряд физических характеристик примесной системы. Взаимодей¬ 
ствие между примесями было рассмотрено в [168]. Все упомянутые 
работы [161—168] ограничены приближением сильно разбавленно¬ 
го сплава. 

Метод когерентного потенциала [160] позволяет рассматривать 
сплав с конечной концентрацией примесей. Можно выделить два 
направления работ, использующих метод когерентного потенциала 
для описания неупорядоченных сплавов. 

Начало первому направлению положила работа [169]. В ней 
была дана теоретическая интерпретация зависимости от концен¬ 
трации средней намагниченности, атомных моментов компонент 
и электронной теплоемкости для сплава ШсЕе^,.. К этому направ¬ 
лению примыкают работы [170—174]. 

Подход Хасегава и Канамори (ХК) основан на использовании 
приближения Хартри — Фока для описания внутриатомной куло¬ 
новской корреляции. В зтом случае гамильтониан (69) записывал¬ 
ся в следующем виде [169]: 


2 Іц а іо а іа +2 Е іа а\ с а іа ,' 

(71) 

Еіа=г | + С7і ( п і-а)- 

(71а) 


Таким образом, неупорядоченность, описываемая в рамках при¬ 
ближения когерентного потенциала, характеризуется двумя пара¬ 
метрами Е Аа и Е в „. Средние числа заполнения (га г _ 0 ) в (71а), 
которые различаются для разных компонент сплава «Пі-о) = п а<у 
или «на, і^А, или В), должно определяться самосогласованным 
образом. Последнее обстоятельство приводит к тому, что не каждая 
элементарная ячейка является электронейтральной и может иметь 
место перенос конечного заряда. 

Для одночастичного гамильтониана (71) применима стандарт¬ 
ная схема метода когерентного потенциала, которую здесь опи¬ 
шем, следуя обозначениям работы [160]. В методе когерентного 
потенциала (СРА) рассматривается одноэлектронный гамильто- 
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ниан следующего вида: 

Я = + = + (72) 

Здесь IV — периодическая часть; Б — сумма случайных вкладов, 
каждый из которых связан с одним узлом. Одноэлектронные свой¬ 
ства сплава вычисляются как средние по ансамблю по всем воз¬ 
можным конфигурациям атомов в решетке. Обычно рассматривают 
усредненную подобным образом одноэлектронную функцию Грина 
О (*): 

(С (г)> = <( 2 - Д- ИГ‘> = (г-\Ѵ- 2)"‘. (73) 

Определим Г-матрицу для данной конфигурации сплава с помо¬ 
щью уравнения 

6 = (0) + (0)Т{0). (74) 

Тогда функциональное уравнение для определения неизвестного 
оператора 2 будет задаваться условием 

(Г [2]> = 0. (75) 

"Уравнение (75) является самосогласованным определением опе¬ 
ратора 2. 

Полагая, что 

П-Е = 2(Дп- 2 п ) = 2 Ѵ п , (76) 

можно ввести локальный оператор рассеяния 

Т п = Ѵ п (1-(С)Ѵ п )~ \ (77) 

С помощью оператора Т п эффективная среда, характеризуемая 
оператором 2, заменяется рассеянием на реальном атоме в данном 
узле п. В методе когерентного потенциала общее условие само- 
согласования (75) заменяется его одноузельным приближением 

(Т п [2]> = 0. (78) 

Таким образом, при этом подходе примесь считается находящейся 
в эффективной среде, функция Грина которой подбирается так, 
чтобы Г-матрица рассеяния на примеси в среднем была равна 
нулю. При этом будем пренебрегать рассеянием парами атомов 
и более крупными кластерами. Метод когерентного потенциала 
точен в атомном пределе, когда перескоки электронов с узла 
на узел очень маловероятны. Сравнение приближений виртуально¬ 
го кристалла, средней Г-матрицы и когерентного потенциала, 
проведенное в [175], показало, что метод когерентного потенциала 
не хуже аппроксимации виртуального кристалла. 
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В методе когерентного потенциала усредненная функция Грина 
неупорядоченной системы {О(Е)) получается из функции Грина 
для идеальной решетки заменой энергии на комплексную величи¬ 
ну. Аналитические свойства величин, вычисляемых в одноузель- 
ном приближении когерентного потенциала, нетривиальны; функ¬ 
ция Грина {О (г)) аналитична всюду, кроме линий разрезов, соот¬ 
ветствующих примесной зоне и зоне основного кристалла. 

Существенно, что в методе когерентного потенциала эффект 
рассеяния электронов вследствие неупорядоченности описывается 
комплексной величиной, а именно когерентным потенциалом. 
С точки зрения квантовой механики в этом нет ничего необычного. 
Напомним, что при многократном рассеянии волны на произволь¬ 
ном ансамбле рассеивателей вводится усредненная по ансамблю 
волновая функция, а потенциал в уравнении Шредингера стано¬ 
вится комплексным [176]. Мнимая часть потенциала описывает 
поглощение вследствие рассеяния. 

Основная характеристика спектра возбуждений системы есть 
плотность состояний на единицу энергии /7 (е). Она определяется 
мнимой частью функции Грина {О (г)) = С СРА . На основе одно¬ 
частичной плотности состояний с помощью метода когерентного 
потенциала можно хорошо описать поведение параметра асферич¬ 
ности у для сплавов Ш, Ре и Со [177]. 

Параметр асферичности является важной характеристикой, 
экспериментально измеряемой с помощью рассеяния медленных 
нейтронов и определяется следующим соотношением: 

'7 = И-е 8 /Г, (79) 

где ц вг; — магнитный момент, определяемый электронами в состоя¬ 
ниях е а -типа; ц — полный спиновый магнитный момент. 

Эксперименты по рассеянию нейтронов показывают, что изме¬ 
ряемые значения у в зависимости от ц очень точно укладываются 
на прямую линию практически для всех сплавов ІЧі, Ре, Со, т. е 

у = а + Ьц. (80) 

Только для чистого № это не выполняется; ущ значительно 
меньше величины, следующей из (80). Возможной причиной 
такого отклонения для чистого № может быть либо влияние 
корреляции электронов, либо специфика одночастичного поведе¬ 
ния системы. В [177] были рассмотрены только одночастичные 
свойства системы в подходе Хасегава и Канамори (71) и показано, 
что для расчета параметра асферичности влияние корреляции 
не очень существенно. Как и в [169], рассматривалась область 
концентраций сплава Ре с при 0^с^0,5. Хасегава и Кана¬ 

мори с помощью метода когерентного потенциала вычислили 
магнитный момент ц и локальные моменты ц (Ш) и ц (Ре). Их 
9-01068 
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результаты хорошо согласуются с экспериментом. Однако надо 
заметить, что они использовали не реальную плотность состояний, 
а сильно идеализированную функцию и проблема решалась с ис¬ 
пользованием многих свободных параметров. 

В [177] впервые была использована реальная теоретическая 
плотность состояний [51, 178] для расчета параметра асферичности 
у. Для точного расчета у необходимо было отдельно учесть е й - 
и ^-состояния. Получить такие раздельные плотности состояний 
весьма сложно из-за сильной гибридизации этих состояний. 
В [177] использовано то обстоятельство, что в точках и на линиях 
высокой симметрии, где гибридизация отсутствует, волновые функ¬ 
ции можно отождествить с е 8 - и ^-состояниями. Предполагалось, 
что количественно поведение волновых функций не сильно изме¬ 
няется при переходе к другим точкам. Используемая теоретиче¬ 
ская плотность состояний состоит из шести подзон, две из них 
связаны с 5 -электронами, а остальные четыре имеют в указанных 
точках и на линиях высокой симметрии поведение плотности 
состояний электронов в І ч - и ^-состояниях. Поэтому можно пред¬ 
положить приближенное разделение плотности состояний на 
составляющие для і 2е - и е 8 -электронов. 

В методе когерентного потенциала выражение для плотности 
состояний в сплаве Ш^Ре,. имеет вид [177] 

= __й. Іт С сра ( 81 ) 


где 


дСРА _ 1 


Лг' 


е—2;—е' ’ 


(82) 


— когерентный потенциал, определяемый из уравнения 

2 г = яД + 2 | (Д-2,)С? РА (е); (83) 

Д описывает сдвиг между атомными уровнями Ре и N1. В [169] 
этот параметр очень сильно зависит от спина (Д*/Д; = 5,6) йот 
концентрации. В [177], напротив, предполагалось, что Д практи¬ 
чески не зависит от этих величин, чтобы последовательно про¬ 
вести учет одночастичных свойств модели. Решение задачи удает¬ 
ся провести без использования свободных параметров. Были 
вычислены плотность состояний В СРА (е) и локальные плотности 
Д 10 ^ 1 (е) и В І0 і' Ѵе (е) для і = і 2е и различных концентраций. 
Полученный на основе этих результатов для Б СРА (е) параметр 
асферичности у показан на рис. 11. Согласие с экспериментом 
хорошее. 

Интересно отметить, что результаты для вычисленных Эльком 
значений ц, ц (Ш) и ц (Ре) оказываются хуже, чем в работе Хасе- 
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гава и Канамори. Возможной причиной этого может быть влия¬ 
ние корреляций на значение р, для описания которой в [169] 
использовали дополнительные свободные параметры. В то же 
время, как видно из рис. 11, поведение параметра асферичности 
хорошо объясняется уже на основе одночастичной плотности 
состояний, оптимально приближенной к реальной. Дальнейшее 
обсуждение подхода Хасегава—-Канамори дано в [179]. 

Другое направление описания неупорядоченных сплавов 
с помощью гамильтониана (69) развивалось в [180, 181]; конкретно 
в [180] рассматривался сплав Р<1 — N1. Подробно проанализиро¬ 
вал различие этих двух подходов Фукуяма [162, 174]. Он показал, 


Рис. 11. Зависимость параметра 
асферичности у от ;концентрации 
х (X — теоретические значения; 

I — эксперимент) [177] 


что в подходе Харриса — Цукермана [180] основное внимание 
сосредоточивается на динамических эффектах кулоновского взаи¬ 
модействия, а пространственным изменением потенциала прене- 
брегается. Поэтому такие одночастичные величины, как локаль¬ 
ная плотность состояний, являются пространственно однородны¬ 
ми, за исключением возможного существования виртуально свя¬ 
занных состояний. Схема является самосогласованной, если имеет 
место равенство е Аа = г Ва в уравнении (69); в этом случае воз¬ 
можно в отличие от (71) учесть некоторые процессы электрон- 
дырочного рассеяния более высокого порядка. 

Различие между подходами Хасегава — Канамори [169, 173, 
179] и Харриса — Цукермана [180] наиболее заметно проявляется 
при рассмотрении коллективных эффектов, в частности, при вычис¬ 
лении спиновой восприимчивости. Это связано с тем, что при по¬ 
строении теории электронных и магнитных свойств неупорядочен¬ 
ных сплавов, описывающихся гамильтонианом (69), необходимо 
учитывать случайное расположение атомов компонент на решет¬ 
ке и влияние кулоновской корреляции электронов на электрон¬ 
ную структуру и физические свойства. Если, как мы видели вы¬ 
ше, одночастичные характеристики сплава (например, параметр 
асферичности у) слабо зависят от корреляционных эффектов, то 
для коллективных свойств правильный учет корреляции более 
существен. 






404 А. Л. КУЗЕМСКИЙ 


9. СПИНОВЫЕ ВОЛНЫ И ИХ УСТОЙЧИВОСТЬ 
В НЕУПОРЯДОЧЕННЫХ ФЕРРОМАГНИТНЫХ МЕТАЛЛИЧЕСКИХ 
СПЛАВАХ 

В последние годы большой интерес вызывает изучение неупру¬ 
гого рассеяния медленных нейтронов в неупорядоченных магнит¬ 
ных металлических сплавах [182—192]. Существует большое число 
экспериментальных данных об изменении спин-волновой жесткос¬ 
ти 77 сплава в зависимости от концентрации. 

Ранние расчеты энергии спиновых волн в сплавах в рамках 
зонной модели основывались на приближении жесткой зоны [52, 
193—195]. В этом подходе сплав описывается как чистый металл 
с соответствующим образом подобранной плотностью электронов 
и зонным расщеплением. По существу, к этому направлению 
примыкает теория Ямада и Шимицу [166]. В 1973 г. Хилл и Эд¬ 
вардс 1196] получили приближенное выражение для спин-волно¬ 
вой жесткости 7? сплава, записанное через одноэлектронные функ¬ 
ции Грина. 

Для модели Хаббарда со случайными параметрами (69) [расчет 
коэффициента спин-волновой жесткости был выполнен Фукуямой 
[172, 174], Ридингер и Научель-Блох [197, 198]. При этом элект- 
рон-электронное взаимодействие учитывалось в приближении 
Хартри — Фока (НР), а разупорядочение — в приближении коге¬ 
рентного потенциала (СРА). Эдвардс и Хилл [199] дали гораздо 
более простой и изящный вывод выражения для спин-волновой 
жесткости в рамках приближения СРА — КРА, используя общие 
выражения (17), (18). Схема расцепления в КРА дана в [200]. 

В [184—192] приведены результаты нейтронных измерений 
спин-волновой жесткости сильно ферромагнитных сплавов на осно¬ 
ве N1 и дана подробная интерпретация результатов с помощью 
зонной модели. В частности, в [191] с помощью неупругого рас¬ 
сеяния нейтронов измерялся коэффициент спиновой жесткости 
сплава Со с для концентраций с = 0,21 и 0,5 (при комнат¬ 
ной температуре) и для концентрации с = 0,05 (при Т = 4,2 
и 293 К). Измерения были проведены с помощью трехосного 
спектрометра. Найдено, что коэффициент жесткости медленно 
уменьшается при увеличении концентрации Со, причем для с = 
= 0,05 величина 77 довольно слабо зависит от температуры. 
Между значениями с = 0,21 и с = 0,5 величина 77 почти не ме¬ 
няется. Теоретическая интерпретация этим результатам сначала 
дается на основе простейшей модели жесткой зоны, а корреляция 
описывается в приближении Хартри — Фока. В этом приближе¬ 
нии рассчитаны О для трех сплавов №Со, №Ге и ШМп. Наилуч¬ 
шее согласие с экспериментом наблюдается для первого сплава, 
наихудшее — для последнего. Это, вероятно, связано с увеличе¬ 
нием различия валентности компонент сплава. Далее используют 
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приближение КРА — СРА. Согласие с экспериментом по сравне¬ 
нию с моделью жесткой зоны улучшается: для сплава ШСо незна¬ 
чительно и более значительно для ШГе и ШМп. Тем не менее, 
как отмечают авторы, приближение КРА — СРА неудовлетвори¬ 
тельно объясняет особенности поведения коэффициента жесткости 
в сплаве Кі х _ с Ге с для с = 0,5. 

Таким образом, дальнейшее улучшение теории должно осно¬ 
вываться на более точном учете межэлектронной корреляции. 
Дело в том, что приближение Хартри — Фока сильно преувеличи¬ 
вает корреляционные эффекты и тенденцию системы к образова¬ 
нию магнитоупорядоченного состояния. Более точный учет корре¬ 
ляции должен ослаблять тенденцию системы к образованию 
магнитного порядка и тем самым уменьшать спин-волновую жест¬ 
кость. Эта идея была реализована Э. Коллей, В. Коллей, А. Л. Ку- 
земским [201]. 

Для учета эффектов корреляции электронов при вычислении 
Б в [202] приведена схема, выходящая за рамки КРА и основан¬ 
ная на когерентном лестничном приближении СБА [203, 204], 
т. е. на самосогласованной комбинации СРА и локального лестнич¬ 
ного приближения [205] в канале частица—частица. Такое Г-мат- 
ричное приближение удобно использовать для сильных коротко¬ 
действующих взаимодействий и малой плотности носителей, так 
что его можно применить к Ш, Рй и РЕ Если зависящую от энер¬ 
гии Г-матрицу заменить эффективным взаимодействием типа 
Канамори [206], как это делалось, например, при вычислении 
парамагнитной восприимчивости [207] и магнитострикции [208], 
то снова приходим к теории спиновых волн в рамках КРА — 
СРА. 

Теория межэлектронной корреляции в локальном лестничном 
приближении в неупорядоченных сплавах, развитая в работах 
[202—204], оказалась очень полезной для исследования спиновых 
волн и их устойчивости как в чистых металлах [209], так и 
в неупорядоченных ферромагнитных металлических сплавах [201, 
210]. В рамках этого подхода вычислим здесь, следуя работе 
[201], величину Б для однозонной модели Хаббарда со случайны¬ 
ми параметрами при нулевой температуре. Корреляционные эф¬ 
фекты определяли с сохранением энергетической зависимости 
Г-матрицы. 

Перепишем гамильтониан Хаббарда со случайными параметра¬ 
ми (69) в следующем виде: 

Н ІѴ) = 2 в к п ка + 2 е і га іа + 2 (84) 

где п ка (п іа ) — оператор чисел заполнения для состояний Блоха 
(Ванье) со спином а; е к — зонная энергия, которая предполагает¬ 
ся независимой от заданной конфигурации {ѵ} сплава. Одно- 
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частичный потенциал и затравочное внутриатомное кулоновское 
взаимодействие Щ принимают случайные значения е ѵ и І7 Ѵ (ѵ = 
'= А, В) в зависимости от того, заполнен ли узел атомом А или 
В. Индекс {ѵ} = {ѵ х , . . ѵ г , . . ѵ^}, где ѵ г - = А, В описы¬ 

вает конфигурацию сплава в целом. 

Энергию спиновых волн со д = Од 2 для кубических кристал¬ 
лов можно определить по полюсу поперечной спиновой восприим¬ 
чивости х + “ (?, со), что приводит к следующему выражению для 
коэффициента жесткости: 

0= ~ 2(4 < ’ 1 >. ’ > ‘ )]- < 85 > 

где 2 ((5\) < ѵ >) с = (п^ — п + ) — намагниченность на узел (п а — 
среднее число электронов на узел); {. . . )( ѵ > — среднее по основ¬ 
ному состоянию при заданной конфигурации {ѵ}; (. . . ) с — сред¬ 
нее по конфигурациям. В работах [196, 199, 200] для сплавов была 
использована формула 

[Пт4г«[5$, д/і„]> {ѵ, > с -1іт Ит хГ (9, со)], 

и д-0 9 ш-0 д-0 ^ 

(86) 

записанная в терминах функции отклика х}~ (ц, со) для спинового 
тока по аналогии с (17). 

Для модели Хаббарда со случайными параметрами (84), фурье- 
образы оператора плотности поперечного спина 5^ (или = 
= ($$) + ) и оператора тока (или = (/^) + ), которые здесь 
не зависят от случайных параметров, имеют вид 

5 ч = 2 а и«к+ д| ; 

Я’?* — у-^г 2 ( е к+Ч~ е к) ак|а к+ 5|, 

где а к а («ко) — оператор рождения (уничтожения) электрона в 
состоянии | ко); N —число узлов решетки. 

Восприимчивости, входящие в (85), (86), могут быть выражены 
[202] с помощью уравнения (87) через причинные функции Грина 
при нулевой температуре следующим образом: 

зГ(ч, ш) =-<«$$, я-_ ч »і ѵ >) с = 
ІІг<ИЛОД<*.*+®; ч)С\"х 
Х(^ + со)А 0 (-д)С| ѵ, (^)}) с ; 



( 88 ) 
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?Хт'(Ч, со) = — <«д/5, д/і ч »і ѵ, ) с = 

= - ■ іг 1 Ш <* г < Л т (Е,Е + со; ц) Сі ѵ) х 

Х(^ + со)Я 1 (- Ч )СГ , (^)}> с , (89) 

где 

Лй}(Е,Я+»; 9) = 

н 

= *«у (ч)-ву } Ц-і^іПи (Я, ^ +со; СО) 2 б|& (I) X 

хЛ^Я, Я + ю; д)Я$(Ё + <о) (а = 0, 1); (90) 

^огу(ч) = ехр( —щК г )б и ; Х іи (ч) = 

= г, ; [ехр( — іцК г ) — ехр (— іцК ; )] ; 
іі) = 1/ЛГ 2 е„ ехр [ік (К* — К,)]. 

Здесь предполагалась локальность только неприводимой вершины 
частица — дырка / и { | } + * (Е, Ё + со; со) == / 4 , (Е, Ё + со; Е + 
+ со, Е ), а шпур означает суммирование (без спина) по одночас¬ 
тичным состояниям. Как видно, формула (86) более удобна для 
вычислений, чем (85), так как уравнение Бете — Солпитера (90) 
для а = 1 можно решить без привлечения дальнейших предполо- 
жений относительно. Разложение ^ и эффективного тока 
Л^> с переворотом спина, входящих в (89) и (90), в первом поряд¬ 
ке по ц и учет кубической симметрии приводит к формулам: 

ХГ (Ч = о, со) = ^ | -|§- <1г {ЛЭД (Е,Е + со) X. 

Х0\ ѵ) (Я + со) і^ ѵ} (Я)» е ; (92) 

Л 1 Л>(Я,Я + со) = .|;,-6у Я + со;со)Х 

X 2 ОД(Я)ЛЙ| П (Я, Ё + со) 6ЙІ (Я + со), (93) 

где введены обозначения }ц= — і^(К г — К^); 

М<7) = Ч‘І; Лін (Я, Я + со; д) = Ч .А[]\(Е,Е + с о). 

Под знаком шпура ] и образуют скалярное произведение. 
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Отделяя диагональные и недиагональные части Л { ^ } в (92) 
и (93), получаем: 

хУ (9 = 0, (о)= 

= ж\^(Ь№Г( Е + <о)№І ѵ итс + ХУ(Ч = О,<0), (94) 

где 

ХГ (Я -= 0, «>) = (2 Л ш (Е, Е + со) X 

хК&Ѵ(Д+о>, Е)У с , (95) 

К/&(Д+«>,Д) = 

= 2 ей?/ (Е + (й) і тп С‘І| ( Е ). (96) 

Поскольку конфигурационное усреднение в (95) выходит за 
рамки СРА, воспользуемся приближением [202] следующего 
типа: (Л 1 і ѵ, -К {ѴІ ) С = (А} Ѵ} ><К {Ѵ> ), так что 

К и (Е + со, Е) = <СІ Ѵ) (Е + со) ІС| Ѵ> (Е)) сіі = 

= 4-2 *МІ (Е + ®) (ВД к е к = 0 (97) 

к 

ХГ (9 = 0, со) = 0, 

т. е. нет вершинных поправок, благодаря симметрии относитель¬ 
но обращения времени. Здесь § ка — когерентная одночастичная 
функция Грина с учетом электрон-электронных корреляций (см. 
ниже). Таким образом, на основе СРА получается формула 

ІХГ (9 = 0, ю) = -4 5 -Ц- 2 $к* (Я+ со) ( Е ) (Ѵкбк) 2 - (98) 

к 

Подставляя (98) и предельное выражение 

1ІШ 4 «[<5д, 9 /-_ Ч ]) {Ѵк )с = 

9-0 9 

= -бж2«^о> <ѵ) >сѴ|е к (99) 

ко 

в 6 (86) и переходя от причинной к запаздывающей г функции 
Грина, получаем 

Д = 6 я(в р в р Іт і <іЕ 42№(Я)-^*(я))*(Ѵквк), (100) 
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где (г — энергия Ферми. Это выражение формально совпадает 
с тем, которое следует из ИРА — СРА [172, 196, 199] на основе 
приближения НР. 

В настоящем расчете, однако, для функции ^ ко используется 
схема СЬА [204]. Тогда корреляционная часть в терминах частич¬ 
но усредненных причинных функций имеет вид 

*тіЛЕ)= { -^СІ-аіЁ) П{Е + Ё), (ѵ = А, В)- (101) 

Гі {Е) = [ту~+ і И (Я) • (102) 

где Т\ — эффект ивна я|двухчастичная вершина. Локальная функ¬ 
ция Грина Вис (г), записанная в виде резольвенты (здесь г — 
комплексная энергия), перенормируется следующим образом: 


Піа ( 2 ) = Р а '(*)/ {Т- [гі (г) - 2 а ( 2 )] Р а ( 2 )}; 

(103) 

е іа (г) = еУ+2шіа(г); 

(104) 

^.(*)=т2^(*)і 

(105) 

Е к 

^ка (г) = (г— е к — 2 а (г)) -1 ; 

(106) 

2 а (г) = се а (г) + (1 — с) Е?(г)- 

- [е $ ( 2 ) - 2 а (г)] Р а (г) [е? (г) - 2 а ( 2 )]; 

(107) 

і йЕ \шР г а (Е). 

(108) 


Здесь 2 а — когерентный потенциал; п — среднее число электро¬ 
нов на узел. В отличие от обычного СРА атомный потенциал 
еіа (г) [в (107) индекс і опущен] приобретает энергетическую 
зависимость благодаря массовому оператору 2 ц і{ (г) за счет кор" 
реляций. В приближении НР решение систем уравнений (103)— 
(108) упрощается, так как вместо (101) и (102) используется 
собственная энергия 2 ц^Ѣ = где п у іа — среднее число 

электронов со спином о на ѵ-узлах, определяемое формулой 

Піо = —— ^ сІЕІтабііа (Е). (109) 

С учетом ^конкретной вершины І и } ѵ /і(Е, Ё + со; со) = 
= — Т { У(Е -\-Е + со) и заменой ѵ ->-{ѵ} находим из уравнений 
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(90), (101), (102) соотношение типа Уорда — Такахаши 
соЛЭДі (Е, Е+і о, ч) 6„-ЛМ„ (Е, Е + со, Ч ) = 

= ехр (- і Ч К г ) СЩ" 1 (Е + ш) - ОЩ~ 1 (Е) ехр( - І Ч К ; ), (110) 

(^ {ѵ,_1 (Е))і]а — (Е еУ) 8 и - і и -Ша (Е) 8 І} . (111) 

Рассмотрим теперь вопрос об устойчивости основного состоя¬ 
ния ферромагнетика относительно спин-волновых возбуждений. 
Для чистых металлов этот вопрос подробно обсужден в [211] 
в зонной модели в приближении НР и в [209] ЬЬА. Условие устой¬ 
чивости основного состояния ферромагнетика 

0 = 0(п*-л*)>О (112) 

можно получить из спектрального представления 

со * 

х-’С, «>- -4- [*>' (ИЗ) 

где спектральная интенсивность /з+ 8 _~ (®) 0 относится к кон¬ 

фигурационно-усредненной системе. Магнонный полюс 

ХроТКч, ©) = К — »+)/(ш — Бд 2 + ів) (114) 

для малых д и со можно выделить из континуума Стонера, посколь¬ 
ку спектральный вес возбуждений пар частица — дырка стремит¬ 
ся к нулю при д ->0. Сравнение (113) и (114) приводит к крите¬ 
рию (112). В данном приближении спин-волновое затухание [вме¬ 
сто 8->-0 в (114)] принимает вид 

^ = Г^ 1т ^ + “ Г (°’ ^ 2 ) = 

= 2 ЬцЗІ* (Р) ІтЗІДц) (Ѵк8 к ) 2 , (115) 

как и в КРА — СРА [174], однако здесь учитывается электрон — 
электрон-рассеяние. 

Рассмотрим теперь численные результаты СЬА по О и сравне¬ 
ние их с экспериментальными данными по рассеянию нейтронов 
для сплавов №Ре и ШРй. Для проведения части расчета Б в ана¬ 
литической форме воспользуемся упрощенными выражениями [212]: 

(116) 

(117) 
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где іѵ — полуширина зоны; ѵ т — величина порядка іѵа\ а — 
параметр решетки. Суммирование по к в (105) приводит с учетом 
невозмущенной плотности состояний (116) к функции Грина 

Р„ (г) = (2/ш) (г а — і V 1 — г!); г 0 = [г — 2 а (г)]/іѵ. (118) 

Чтобы (118) было однозначным, выбираем ту ветвь в плоскости 
2 „ с разрезом вдоль действительной оси от —1 до +1, где квадрат¬ 
ный корень положителен на верхнем берегу разреза. 
Переписывая (100) в виде 

В = бя^ 1 -^) Іт 1 йЕ [П Н V' Е > + П П Е ) ~ 2П Н (Я, Е )], 

(119) 

где введены обозначения П рр < (Е, Е) — П оа !(і? + , Е + ),\Е + = Е-\- Ю и 
Паа» (*, *') = ■4- 2 3*0 (2) &ко я (О (Ѵ*е к ) 2 , (120) 


затем применяя теорему о вычетах при суммировании по к 
в (120) с учетом аппроксимации (117), найдем 


П-=^(з г * — 4 — ЗЦ.УТ34); 2р = Д+ -^ . ^ ; 

гтгг 2а= т ( 3 , . (1-^|) 3/2 -(1-# /2 ^ 

п н = ѵ (4+4+ 2 і г і — у+і --) • 


Следовательно, 


д - ^(.І— ,) 1т [«[(ч—*>•- 
-іП^Д(з*,+ 2 1рУ" )- 
іУРЛПа..- 




( 121 ) 

( 122 ) 


(123) 

Скалярная статическая электропроводность о вычисляется 
в том же приближении, что и И. В результате придем к модифици- 
формуле Кубо — Гринвуда 
„ ет 


' 6яУ 


2 [Щ“(р)-ВеП"(ц)]=2°- 


(124) 


которая содержит функции Грина, перенормированные за счет 
ѳлектронныхі корреляций в рамках СЬА. Здесь П рр (р) = 
== П а0 (р + , р~), р“ = р — Ю; V — объем системы; е — единич- 


г 4 



Рис. 12. Зависимость коэффициента спинволновой жесткости И < 0 (неста¬ 
бильный случай) (а); статической электропроводности о а , а (б); плотности 
электронов п а (в); намагниченности тУ, т (г) от ІІ Ѵ (ѵ =А,В ); приближение 
Хартри — Фока; параметры (га, е А , е в , с, п) = (1, —0,8; 0; 0,4; 0,4). 



Рис. 13. Зависимость парциальной намагниченности тУ (о); коэффициента 
жесткости В (б), энергии Ферми ц (в), электропроводности о а , а (г) от п; 
приближение Хартри — Фока; параметры (га, е А , г в , Ь ,А , Ѵ в ) = (1; 0,2; 
0; 4; 3) (а, б, в): 

сплошные линии соответствуют составу с = 0,1; штриховые ■» с = 0,3 
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ный заряд. Вводя (121) и (122) с заменой ЩЗ = П [г, 

—кг а , 2 + в (124), получаем зависящую от спина электро¬ 

проводность (при Іт (ц + ) < 0): 


X {і Ѵі - 1% (і (1 - 2|) + -і- 2 0 Іт Ео (р + )) } ] , (125) 


где 

2 а =[ц+- 2 а (ц+)]/ц;; о = е 2 ^ІѴ/( Зл 2 Ѵ). (126) 


Численный анализ проводится следующим образом: сначала 
выбираются параметры ю, е А , е в = 0, Ѵ А , Ѵ в , с, п, затем решает¬ 
ся самосогласованная система уравнений (101)—(108) с функцией 



Рис. 14. Зависимость спинволновой жесткости О в РеоШ^е от состава. Зна¬ 
чения В ( х ), рассчитанные в когерентном лестничном приближении при 
к е А ,е в , Ѵ А , 17 в , п) = (0,5; —0,24; 0; 2,66; 3,4; 0,6) [(<?) абс. ед.— 2 іѵ = 
= 4,15 эВ, а = 4 А]: 

а — расчет: 1 — [196], 2 — [198], з — [199], 4 — [200] б — эксперимент: □ — [99], 

Д — [184], О — [185] 


Грина (118). Полученные результаты используются для нахож¬ 
дения Ь по (123). 

Переходная область между слабым и сильным ферромагнетиз¬ 
мом в зависимости от внутриатомного отталкивания 1/ ѵ , рас¬ 
сматриваемого здесь в приближении НГ, показана на рис. 12. 
В частности, в [213] были исследованы решения с параметрами 
Ѵ А — 2, Ѵ в = 0 и Ѵ А = 2, Ѵ в = 2, но спиновые волны не рас¬ 
сматривались. Расчеты коэффициента жесткости В (в единицах 
й 0 = (1/9) ы?а г ) в ВРА — СРА (см. рис. 12, а) указывают на неста¬ 
бильность ферромагнитного состояния относительно возбуждения 
спиновых волн. Парциальные л полная средние намагниченности 
тУ = п\ — п\ жпг = тц — изображены на рис. 12, г. Ограни¬ 
ченная область стабильного (в основном, насыщенного) ферромаг¬ 
нетизма Б > 0 дана на рис. 13, б, ш ѵ >0 (см. рис. 13, а), рас- 
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считанная в хартри-фоковском приближении в зависимости от 
плотности электронов п. Нуль Б при наименьшем значении 
п соответствует приблизительно критерию типа Стонера [204], 
тогда как другой нуль говорит об изменении типа магнитного 
порядка. 

Как видно из рис. 14, большое количество нейтронных данных 
для О имеется для сплавов ШГе. В качестве проверки были учтены 



Рис. 15. Зависимость спин-волновой жесткости В ( х) (а); Б в абс. ед. (б), пар¬ 
циальной намагниченности тУ (в); эффективные кулоновские взаимодействия 
Г ѵ , Г (г) от с в Р<1 с рассчитанные в когерентном лестничном прибли¬ 

жении 

(2і») р<1 , (2») МІ , Е РЙ — е Р№ , П* 1 , У 141 =(6,05; 4,15; 0,3; 9,17; 14,11) эВ, п РЙ = 0,4; 
„N1 = о,6; а — —’-расчет [214]; • — эксперимент (ср. [214]) 

межэлектронные корреляции, хотя данная схема более удобна 
для описания компоненты N1 с 0,6 дырками на атом в й-зоне, чем 
компоненты Ре с высокой дырочной плотностью (см., например, 
на малое значение Б при с = 0,4 на рис. 14, б). Сравнение расче¬ 
тов Б (с) на основе СЬА с результатами, полученными в рамках 
КРА — СРА [196, 198], теории «жесткой зоны» [199], подхода 
[200], а также с данными по неупругому рассеянию нейтронов 
[199, 184, 185] дано на рис. 14. Коэффициент жесткости Б для 
чистого Ш оказывается близким к Рмі= 555 МэВ-А 2 , измеренной 
при 4,2 К [99]. 
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Расчеты для сплавов №Рй были проведены в СЬА и эти резуль¬ 
таты (рис. 15) сравниваются с данными работы [214]. Параметры 
для чистых систем выбирались на основе [208]. Сплав форми¬ 
руется сп — сп ра + (1 — с) п х1 , а различные интегралы пере¬ 
скока учитываются через ширину зоны 2 и> — с (2іѵ) рі -(- (1 — 
— с) (2и?) Кі . Заметим, что в приведенных единицах и рл , 17 т , 
е Рй — е КІ для всех с приводятся к шкале 2и> = 1. Затравочные 
ІІ Ѵ перенормируются самосогласованным образом, в результате 
имеем двухчастичные вершины Т'І (Е Е). В частности, Г ѵ = 
= Ті (2|х) и Г = сГ А 4- (1 — с) Г в в приведенных единицах 
изображены на рис. 15, г. 

Таким образом, численные результаты, полученные в рамках 
СЬА, указывают на влияние электрон-электронных корреляций 
на энергию магнонов в длинноволновом пределе = Бд*. Не¬ 
смотря на однозонное приближение в модели Хаббарда с упро¬ 
щенной структурой зоны и диагональным беспорядком, для спла¬ 
вов на основе Ш получены физически разумные значения Б. Для 
зависящих от энергии двухчастичных вершин предполагается 
локальность, что позволяет сохранить одноузельный характер 
СРА. Спин-волновое затухание оказывается малым, порядка 
< 7 4 . В этом смысле предложенный в [201, 202, 209] самосогласо¬ 
ванный метод для нахождения стабильного ферромагнетизма 
в сплавах переходных металлов представляется интересным и пер¬ 
спективным. Данный метод может также оказаться полезным 
при исследовании изменения спин-волновой жесткости сплавов 
ферромагнитных переходных металлов в зависимости от степени 
упорядочения сплава. 


ЗАКЛЮЧЕНИЕ 

Исследование упругого и неупругого рассеяния нейтронов 
в магнитных переходных металлах и их сплавах по-прежнему оста¬ 
ется одним из наиболее интересных с принципиальной точки зре¬ 
ния направлений магнитной нейтронографии. Дальнейший про¬ 
гресс в теоретическом описании магнитных свойств переходных 
металлов и их сплавов тесно связан с экспериментальными иссле¬ 
дованиями, которые должны дать точные и надежные сведения 
о распределении зарядовой и спиновой плотностей, спектре 
магнитных возбуждений, обменных связях, силах межатомного 
сцепления и других физических характеристиках кристаллов. 
Возможности магнитной нейтронографии предоставляют уникаль¬ 
ные экспериментальные данные, в частности, о спектре магнитных 
возбуждений. 

В настоящей работе показано, что широкий круг явлений 
в ферромагнитных переходных металлах и сплавах, изучаемых 
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с помощью рассеяния медленных нейтронов, можно последова¬ 
тельно описать с единой точки зрения, используя модель Хаббар¬ 
да. Спектр магнитных возбуждений модели Хаббарда отражает 
двойственный характер поведения магнитоактивных й-электронов 
в ферромагнитных переходных металлах. Исследование свойств 
ротационной инвариантности гамильтониана Хаббарда показы¬ 
вает, что наличие спин-волнового акустического полюса у маг¬ 
нитной восприимчивости является прямым следствием ротацион¬ 
ной симметрии системы, как и для модели Гейзенберга. Таким 
образом, акустическая спин-волновая ветвь — отражение опре¬ 
деленной степени локализации ^-электронов; характерная вели¬ 
чина 2), определяющая спин-волновую жесткость, непосредствен¬ 
но измеряется в нейтронных экспериментах. В отличие от зонной 
модели ферромагнетизма, основанной на предположении о полной 
коллективизации магнитных электронов, модель Хаббарда позво¬ 
ляет учитывать корреляционные эффекты. Поэтому с помощью 
модели Хаббарда можно в принципе описать (при достаточно 
точном учете корреляционных эффектов) большую совокупность 
не только магнитных, но и электрических свойств кристаллов. 
Уже простое приближение хаотических фаз позволяет достаточно 
удовлетворительно рассчитывать спектр магнитных возбуждений 
модели Хаббарда; данный спектр, помимо акустической спин- 
волновой ветви, содержит континуум одночастичных стонеровских 
возбуждений. Наличие стонеровского континуума — проявление 
делокализации магнитных электронов. 

Поскольку стонеровские возбуждения не возникают в модели 
Гейзенберга, их обнаружение и детальное исследование с по¬ 
мощью метода рассеяния нейтронов составляет одну из актуаль¬ 
нейших проблем физики магнитного состояния. Более того, сле¬ 
дует ожидать, что измерения в «оптическом» диапазоне дадут более 
интересные результаты, чем те, которые можно предполагать на ос¬ 
новании однозонной модели Хаббарда. Последняя модель, отра¬ 
жая основные черты поведения системы, сильно упрощена. Мо¬ 
дель Хаббарда допускает естественные обобщения, позволяющие 
учитывать орбитальное вырождение й-уровней, (а — ^-гибриди¬ 
зацию, электрон-фононное и электрон-магнонное взаимодействие 
и т. п. Мы показали, что уже простейшие обобщения модели Хаб¬ 
барда дают значительно более {богатый спектр, причем именно 
в «оптическом» диапазоне. Спин-волновой акустический’полюс со¬ 
храняется, если сохраняется ротационная инвариантность моди¬ 
фицированного гамильтониана. “і 

Модель Хаббарда чрезвычайно эффективна для теоретического 
описания электрических и магнитных свойств сплавов переход¬ 
ных металлов. Введение случайных параметров позволяет смоде¬ 
лировать "неупорядоченный сплав ферромагнитных переходных 
металлов. Модель Хаббарда со случайными параметрами содер- 
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жит в качестве предельных случаев большинство известных моде¬ 
лей систем с примесями. Использование приближения когерентно¬ 
го потенциала позволяет рассчитывать целый ряд физических 
величин, наблюдаемых в нейтронных экспериментах, в широком 
интервале концентраций примеси. При расчете коллективных 
характеристик сплава, таких, как спектр спин-волновых возбуж¬ 
дений и статическая электропроводность, важно последовательно 
учитывать корреляционные эффекты. В данной работе было пока¬ 
зано, что энергию длинноволновых спиновых волн в сплавах 
ферромагнитных переходных металлов можно вычислить, исполь¬ 
зуя приближение когерентного потенциала для описания неупоря¬ 
доченности и общие соображения о ротационной инвариантности 
системы. Метод расчета, развитый Э. Коллей и В. Коллей, исполь¬ 
зует идеи микроскопической теории ферми-жидкости. При нуле¬ 
вой температуре в когерентном горизонтальном лестничном при¬ 
ближении для модели Хаббарда со случайными параметрами про¬ 
ведено самосогласованное вычисление ренормировки коэффициен¬ 
та спин-волновой жесткости за счет электрон-электронных кор¬ 
реляций. Коэффициент П получен численным образом и исполь¬ 
зуется для определения устойчивости ферромагнитного состояния. 
Сравнение полученных результатов с зкспериментальными дан¬ 
ными по рассеянию нейтронов для сплавов на основе № показы¬ 
вает, что в области применимости наши расчеты хорошо согласуют¬ 
ся с экспериментом. При низкой плотности предлагаемый подход 
позволяет более точно учитывать корреляционные эффекты, чем 
в приближении хаотических фаз. 

Теория магнитных явлений в ферромагнитных переходных ме¬ 
таллах и их сплавах ни в коей мере не завершена. Теория рас¬ 
сеяния медленных нейтронов в ферромагнитных переходных ме¬ 
таллах, основанная на приближении хаотических фаз, нуждается 
в ряде уточнений [136, 215]. Спектр стонеровских возбуждений 
Йсо ч = 8 к+д — е к + А в более тойном приближении [216] имеет 
вид 

/гсо ч = е к+ч — е к + 5 (к + ц). 

Насколько хорошо применимо приближение хаотических фаз, 
следует судить по результатам эксперимента. В настоящее время 
также есть веские основания полагать, что более детальная теория 
неупругого рассеяния нейтронов в ферромагнитных переходных 
металлах должна учитывать зависимость кулоновской корреляции 
от квазиимпульса и орбитальное вырождение. 

Большой интерес представляет дальнейшее теоретическое и 
экспериментальное исследование ферромагнитных сплавов пере¬ 
ходных металлов. В последнее время рассеяние медленных нейтро¬ 
нов применяется для измерения коэффициента спин-волновой 
жесткости в ферромагнитных частично упорядоченных сплавах 
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[185, 217, 218]. В недавних теоретических работах по исследова¬ 
нию частично упорядоченных ферромагнитных сплавов [219—224] 
используются многие идеи, описанные в настоящем обзоре, в част¬ 
ности в [222], которая является дальнейшим развитием работы 
[199], разработана теория, позволяющая вычислять изменения 
коэффициента спин-волновой жесткости в зависимости от степени 
атомного упорядочения в ферромагнитных сплавах переходных 
металлов типа А„ і5 В 0і5 . 

Предполагалось, что решетку сплава можно представить в ви¬ 
де двух взаимопроникающих подрешеток компонент А и В, так 
что все ближайшие соседи атомов из подрешеток А принадлежат 
подрешетке В, и наоборот. Степень атомного упорядочения харак¬ 
теризуется величиной р, которая представляет собой вероятность 
того, что атомы А (В) принадлежат подрешетке А (В). Вероят¬ 
ность того, что атомы А (В) принадлежат подрешетке В (А), рав¬ 
няется (1 — р). Упорядоченное и полностью разупорядоченное 
состояния характеризуются значениями р = 1ир = 0,5 соответ¬ 
ственно. Параметр р связан с обычно используемым параметром 
дальнего порядка 5 уравнением 8 = 2 р — 1. Электронная под¬ 
система описывается гамильтонианом Хаббарда (69). Для описа¬ 
ния структурной неупорядоченности использовался метод коге¬ 
рентного потенциала. Стандартное одноузельное рассмотрение 
позволяет получить систему уравнений, определяющую коге¬ 
рентный потенциал, числа заполнения и эффективный потенциал. 
Для вычисления энергии спиновых волн, как и в работе [199], 
использовались общие формулы (17), (86). Межэлектронная кор¬ 
реляция учитывалась в приближении хаотических фаз. Конфигу¬ 
рационное усреднение проводилось с учетом подхода Белицкого 
[ 212 ]. 

Для простой модели сплава проведены численные расчеты 
коэффициента спин-волновой жесткости в зависимости от пара¬ 
метра атомного порядка р при нулевой температуре. Выбор пара¬ 
метров соответствует случаю сплава ШРе, хотя плотность состоя¬ 
ний используется сл и шком упрощенная. Последнее обстоятельство 
затрудняет прямое сравнение развитой теории с экспериментом, 
тем более что нейтронные измерения зависимости Б от степени 
атомного порядка проводились для № 3 Ре и РД 3 Ре. Тем не менее 
вычисленное поведение Б в целом отвечает эксперименту, т. е. 
увеличение степени неупорядоченности приводит к спин-волно- 
вому «размягчению» системы. Дальнейшее теоретическое и экспе¬ 
риментальное изучение подобных систем весьма интересно. Ме¬ 
тод неупругого рассеяния нейтронов в данной задаче — единст¬ 
венный метод проверки теоретических результатов. 

Дальнейшее всестороннее теоретическое и экспериментальное 
изучение ферромагнитных переходных металлов и их сплавов 
приведет к углублению наших представлений о природе магне- 
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тизма, уточнению модельных теоретических представлений и при¬ 
менимости приближенных данных о их магнитных характеристи¬ 
ках, что имеет важное значение для приложений и понимания 
магнитных и электронных свойств этой обширной группы веществ. 

В заключение мне хотелось бы выразить благодарность моим 
соавторам Л. Черу, Э. Коллей, В. Коллей за интересные и стиму¬ 
лирующие обсуждения. Я искренне благодарен моим товарищам 
по работе Н. М. Планиде, В. Л. Аксенову и В. Б. Приезжеву за 
их большую помощь и внимание. 
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